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Abstract. We consider the spreading of a thin two-dimensional droplet on a solid sub- 
strate. We use a model for viscous fluids where the evolution is governed by Darcy's Law. 
At the triple point where air and liquid meet the solid substrate, the Uquid assumes a con- 
stant, non-zero contact angle (partial wetting). We show local and global well-posedness 
of this free boundary problem in the presence of the moving contact point. Our estimates 
are uniform in the contact angle assumed by the liquid at the contact point. In the so-called 
lubrication approximation (long-wave limit) we show that the solutions converge to the 
solution of a one-dimensional degenerate parabolic fourth order equation which belongs to 
a family of thin-film equations. The main technical difficulty is to describe the evolution 
of the non-smooth domain and to identify suitable spaces that capture the transition to the 
asymptotic model uniformly in the small parameter e. 
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1. Introduction and model 

In the past years, the theory of fluid systems in the presence of a free boundary has 
been developed by many important works. Usually, in these problems, the interface (or 
free boundary) separates two phases of the fluid system. Among the large literature, such 
work has been addressed e.g. in [31, 30, 13, 39, 44, 12] for local existence results, [42, 19] 
for global existence results, [10] for the study of blow-up, [2, 32] for asymptotic limits. In 
this paper, we are interested in the situation of a fluid evolution in the presence of a contact 
point where three phases meet, namely a contact point between air, Uquid and soUd, see 
Fig. 1. One example is the flow in a Hele-Shaw cell, where the liquid touches the lateral 
boundary of the Hele-Shaw cell. 

The Hele-Shaw model describes the evolution of a hquid between the plates of a Hele- 
Shaw cell. In general, the surface tension-driven Hele-Shaw flow is given by 

'Ap = in Q.{t), 

, p = YK on dQ.{t), (1.1) 

V = Vp-n on dQ.{t), 
1 
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where the evolving domain Q.{t) C describes the region occupied by the fluid. The 
velocity of the fluid is described by Darcy's Law U = —Vp. In particular the normal 
velocity V of the fluid interface is described by V = Wp-n. The parameter 7 describes 
the surface tension between air and liquid. Next to its interpretation as the flow in a Hele- 
Shaw cell, the fluid evolutions govemed by Darcy's Law appear in a wide range of physical 
models. One example is the flow of a liquid through a porous mediimi, see [6]. Other 
situations which can be modeled by (1.1) are crystal growth or dissolution, directional 
solidification or melting, electrochemical machining or forming [41, 38, 33]. In the last two 
decades, well-posedness of ( 1 . 1 ) has been investigated: Short- time existence and regularity 
of solutions of (1.1) have been proved in [15, 24, 17, 18] andProkert [35]. Global existence 
for initial data close to the sphere has been shown in [11]. The case of zero surface tension, 
7 = has been considered e.g. in [40, 1]. 

Clearly, the normal component of the velocity is zero at the hquid-sohd interface. We 
assume that the Hele-Shaw cell is described by the half-space H = Rx (0,°°). At the point 
where air, hquid and sohd meet, we assume that the liquid assumes a static (microscopic) 
contact angle. This contact angle 6 is determined by Young's Law [43], i.e. 7cos0 = 
YsG ^ YsL where the parameter describe the surface tensions between the three phases: 7 
(air, liquid), ysL (solid, liquid) and ysG (air, solid). This leads to the following model: 



■Ap = 
p = JK 

Py = 
1.7COS0 = JSG-JSL 



mQ.{t), 
ondQ.{t)nH, 
ondQ.{t)ndH, 
on d{dQ.{t)r\dH), 



(1.2) 



see Fig. 1. The evolution then can also be interpreted as the spreading of a droplet on a 
plate. A well-posedness result for (1.2) in Holder spaces has been given by Bazalyi and 
Friedman in [5, 4]. However, in their analysis the conditions on the initial data are too 
restrictive to allow for movement of the triple point (and thus for spreading of the droplet). 
Our first main result is a well-posedness result for this free boundary problem in a much 
wider class of weighted Sobolev spaces. In particular, our result seems to be the first result 
which allows for movement of the triple point. 

The second aim of this work is to show the convergence of solutions to a reduced model 
in the so called lubrication approximation regime or long wave approximation. More pre- 
cisely, let us assume that typical vertical length scales are of order e while horizontal length 
scales are of order 1. In particular, the angle assumed at the contact point is of order e: 

^ ^^^l-cos0= ^^±^^^. (1.3) 
2 2 7 
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The limit model is a special form of the thin-film equation. Assuming that the height of 
the droplet is described by the graph h{t,x), the evolution is given by 

'ht + Y{hh^)^ = m{h>0}, 

h = 0, \hx\=e, ond{h>0}, (1.4) 

, V = Yh^ax on d{h > 0} 

and where V is the velocity of the moving contact points d{h>Q}. Formal derivations of 
lubrication models of type (1.4) have been e.g. given in [37]. We prove convergence of so- 
lutions of (1.2) to solutions of (1.4). This is the first rigorous lubrication approximation in 
the case of partial wetting (non-zero contact angle). Furthermore, it is the first convergence 
result in the framework of classical solutions. A rigorous lubrication approximation in the 
framework of weak solutions has been done by GiacomelU and Otto [23]. Their approach 
is quite different to ours: In particular, their result does not include well-posedness for the 
initial model. Instead, the authors prove convergence to the limit model by only minimal 
energy bounds, assuming existence of smooth solutions. In particular, the bounds derived 
in [23] do not capture the slope of the profile at the contact point. Indeed, the techniques 
used in [23] do not seem to be apphcable for the case of a non-zero contact angle at the 
moving contact line as considered in this work. 

The main work lies in the derivation of bounds which are uniform in the parameter 
£ > 0. We give a short sketch of the strategy of our proof: In order to perform the transition 
from (1.2) to (1.4), we first express (1.2) as a nonlocal evolution problem in terms of the 
profile function h{t,x). The corresponding equation can be seen as a nonlocal parabolic 
evolution problem of third order. Equation (1.4) on the other hand is a local fourth order 
degenerate parabolic equation. As the considered models are higher order equations, the 
maximum principle cannot be used. Instead we rely on their dissipative structures. Indeed, 
solutions of (1.2) satisfy 

^((\dClnH\)-a\dQ.ndH)) = - ^ [f I'^pfdxdy, (1.5) 
at YJJ^{t) 

where a = cos(arctane), see e.g. [22]. The dissipation relation for solutions of (1.4) is 

^(/ hldx + e^\{h>0}\) = -y[ hhl^dx. (1.6) 
at JwL Jm. 

One of the core issues of the analysis is to find suitable norms which allow for uniform 
bounds in the limit e — >^ 0. In [26], we have investigated the linearizations of (1.2) and 
(1.4). The analysis in this work suggests to use sums of weighted Sobolev norms of the 
type 

" }=}++}- 

where f = ~ I and where k,5 >0. In the Umit e — >■ 0, this norm turns from a sum 
of weighted Sobolev norms of order 4^ -|- 5 and 3^-1-5 to a weighted Sobolev norm of 
order 4k + 6 (first term on the right hand side of (1.7)). This transition in the character of 
the norm is reflected by a transition in the character of the equation: In the limit e — )• 0, 
the model (understood as an evolution equation for the profile fimction, cf. (2.16), (2.19)) 
changes 

• from a third order to a fourth order evolution equation and 

• from a non-degenerate parabolic to a degenerate parabolic equation. 
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In particular, we will use norms of type (1.7) with 5 = 1. This seems to be the smallest 
integer value that is sufficient to control the nonlinearity of the problem (the norms (1.7) 
are stronger for larger 5 as follows from Hardy's inequaUty). We also need to choose 
norms which control the pressure p. As we will see, the norms for the pressure do not have 
a real space representation, but are rather described in terms of the Mellin transformed 
function. In fact, the choice of suitable norms for the pressure turns out to be delicate 
in order to obtain uniform bounds in the small e parameter. We use radial variables with 
respect to the moving contact point at the origin of the coordinate system. The norms, we 
use are weighted Sobolev norms of supremum-type in the angular direction (in real space 
variables) and of L^-type in radial direction (in frequency variables). 

Structure of the paper: In Section 2, we transform the problem on a fixed domain and 
we define the norms to control the profile and the pressure. The main results of this work 
are stated and discussed in Section 3. In Section 4, we give an overview for the proofs of 
the main theorems. In Section 5, we prove estimates for weighted spaces. In Section 6, we 
derive estimates for the nonlinear operator for a droplet supported in half-space. In Section 
7, we derive corresponding localized estimates for compactly supported droplets. 

2. Setting and norms 

By a change of dependent and independent coordinates, we reformulate the problem on 
a fixed domain. We then formulate (1.2) as a nonlocal evolution equation in terms of the 
profile function h. We also introduce norms to control profile and pressure. 

2.1. Transformation onto a fixed domain. We define the nonlinear operator of Dirichlet- 
Neumann type by 

r Ap = inn(/), 

B^K{x) = ^\+hl {dnp\y^h{t,x)^ where \p = 7K ond^£l{t), (2.1) 

\ Py=Q 3XdoQ.{t), 

where d\Sl{t) := dQ.{t) n {>■ > 0} and dQQ.{t) := dQ.{t) n {y = 0} With the assumption 
that the free boundary moves with the velocity, we have ht =w — vhx, where u = (v, w) 
is the velocity of the liquid. Suppose that the support of the droplet stays an interval for 
some time, i.e. supp/!(f) = (i_(f),s+(f)) for t G (0,t). The evolution (1.2) can then be 
equivalently written as nonlocal evolution for the profile h by 

ht+B''K^O, forxe (i_(f),i-(f)) (2-2) 

with boundary conditions h\^x=s± = 0' |/j;c|x=s±l = £ ^iid s±{t) = — i^x\x=s±- By (2.1) and 
since \hx\ = e at the triple point, the movement of the triple point is given by 

dts± = -Px\x=s^ = -rK-x|x=.v± = 7((1 + e^)"^ hxxx - 3e(l + e^)-^ hi) . (2.3) 

Indeed, (2.3) follows since p{t,x,h{t,x)) = K{t,x) and hence Px + hxPy = fCc at x = s{t) 
and we conclude since py = Ofory = 0. We will assume that the support of h at initial time 
is given by (0, 1). We next rescale time and space to get 0(1) quantities; furthermore we 
fix the position of the moving contact point s{t) by using moving coordinates: We set 

M(f) = D{t) = (2.4) 

We introduce the new variables {x,y,t} by 

x-M(t) 1 - J . f' ^ J .^.x 
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cf. [3, 20]. In particular, dxX — g, dxi — 0, = ^ and 

where the dot denotes differentiation in i. The dependent quantities h and p are defined by 
/j(?,x) =eD(?)/i(f,jc), D{f)p{t,x,i) = e'Ypi^,x,f). (2.6) 

The transformed evolution — after multiplication by /ey — is given by 

{D^.-{{\-x)s-+xs+)h-Bl{{l+e^hl)-'^hxx) = Q forxe(0,l) (2.7) 
with boundary conditions h = Q and = 1 at x = 0, 1. The operator Sg is given by 

{Ag/? = in n, 
p = n ondiQ, (2.8) 

Py = on 

and where Ag = + (1)^9^. Here, Q. = {{{xj) : x e (0, l),y e iO,h{t,x)))}. We will 
also use the notation f = d\Q. for the air-liquid interface and Fq = doQ. = (0, 1) for the 
liquid-sohd interface. By Proposition 6.7, is well-defined. In the following, we skip 
the tilde's in our notation. Let h* = x(l — x) and let f* = h* = l — 2x. Note that h* is an 
approximation of the stationary solution for the Darcy flow which is the half-circle. We 
set / = {h — h*)x = hx — f* . We also use the notation b{ := B'^. This yields the following 
evolution model, defined on the fixed domain := (0, t) x (0, 1), 

' Csf := (Df), - [((1 -x)i_ +xi+) (/ + /*) +Bi( ^:'^/\. 

L ^ (l + e2(/- 2)2)2 ■ 

/ = atx = 0,l, (2.9) 

I 

fdx = Q. 



0, 







The contact point positions s± are given by S- (0) = 0, 5+ (0) = 1 together with the ODE 

5±(f) f(l+e2)-i/„-3e3(l+e2)-V2(/^ + 2)2)^ ; (2.10) 



lx=oa 

D[t) is accordingly defined by (2.4). We have transformed the equation onto a fixed domain 
at the cost of the non-local operator b{ as and the nonlocal terms s±, D. The analogous 
transformations for the thin-film equation (1.4) yield 

Cof := {Df), - [((1 -x)i_ +xs+){f + f) + ((£/a?x)/«,) = 0, 
/ = 0, atx = 0,l, (2.11) 

fdx = 0. 

For e = 0, the functions s± are defined by s_ (0) = and s+ (0) = 1 together with the ODE 
s±{t) = fxx\x=o,u D{t) is defined accordingly by (2.4). 

Due to the degeneracy of the evolution equation at the free boundary, special attention 
needs to be directed at the boundary conditions: The above transformations are such that 
the boundary conditions h = 0,\hx\=e respectively, are equivalent to the integral/boundary 
conditions // = 0, / = respectively for the transformed function /. In the subsequent 
part of this work, we construct / by a Lax-Milgram argument such that / = indeed holds 
at the boundary. The integral condition / fdx = holds automatically for sufficiently 
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smooth solutions. Indeed, suppose that / fdx = is satisfied initially. For the Darcy flow, 
with the analogous calculation as for (2.3), it foUows that 



d 
dt 



A 

^{l+e^{f+f*)2)i 

-S++S- 0. 



Ho/dx = -\(il-x)s-+xs+)if + f*)-Bi( ^] 



Therefore, / D fdx = JD fmdx = for all times (as long as D 7^ 0) and hence also / D fdx = 
f Dfi^dx = 0. An analogous calculation applies also for the thin-film equation. Note that 
instead of proposing a fixed contact angle which implies the speed of propagation, derived 
in (2.3)), another option would be to consider a model with dynamic contact angle while 
imposing a law relating contact angle and speed of propagation, see e.g. [36]. 

The case of an infinite wedge. Near the moving contact lines, the region occupied by 
the liquid approximately has the shape of a wedge. This motivates to linearize the evolution 
equation around an infinite wedge. We hence assume that h{t,x) « e(x — s(f )) for .s-(f ) G R. 
We describe how the problem is transformed onto the wedge, a more detailed derivation is 
given in [26]. Analogously to (2.5), the new variables are defined by 

x — s{t) = X, y = ey, t = *- and p = sp, h = eh. (2.12) 

Correspondingly to (2.7), we get 

k-sh + Bl{{l+eW^)-hja) = forxe(0,-). (2.13) 

in the following, we omit the in the notation. Here, s{t) is defined as in (2.3). We set 
f :=hx — 1. Taking one spatial derivative, of (2.13), we get 

f-[sf+Bl{f,{l+e\\+ff)-l)]^ = (2.14) 

with the single boundary condition f\x=o = 0. We introduce the notation = {{x,y) :x> 
Q,0 <y <h},Yf = {{x,y) : x > 0,y = h}. The linear operator Bg := Bg is given by 

A^q = mK, 

HyJiy=,, ^ 



Be'n{t,x) = {-qx + ilf qy),^^^, where { q = r] ondiK, (2.15) 



qy = on ^0^) 

where K := and F := F^. Equation (2.14) can then be equivalently expressed as 

■f+Aef = Neif) forxe(0,oo), 
^ A f n (2.16) 

/ = for X = 0. 

The main (linear) part of (2.14) is given by the operator 

A, := ~ (l+e2)-2 d,B,d,. (2.17) 
The remaining terms in (2.14) are combined in the nonlinear operator Ne{f) — Ne{f,f), 

The first term on the right hand side of (2. 1 8) is related to the movement of the triple point. 
The second and third term sum describe the error which appears by replacing the domain 
Rf by K and by replacing the curvature with fx. Analogously, for the thin-film equation 
we apply the coordinate transform x = x — s{t). In the new coordinates, we obtain 

(ft+Aof = No{f) forxe(0,~), 
I / = forjce (0,00). 
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The linear and nonlinear part of the equation are given by 

Aof ■= {xf„),„ No{(pJ) := -{Io(pdx{f„-f,,i,=o))^ (2.20) 

and No{f) = No{f,f). We also write Aq = —dxBodx and Bq = —dxxdx. Observe that 
No{(p,f) is bilinear, while Ne{(p,f) is neither linear in the first nor in the second argument. 
Also notice that (2.19) has the scaling invariance (x,f ,/) i— >■ {Xx,X^t,f ). 

2.2. Norms for the profile. The initial problem (1.2) is non-degenerate parabolic on a 
non-smooth moving domain, the Umit problem (1.4) is degenerate parabolic on a smooth 
domain. We use weighted Sobolev type spaces to capture the transition between these two 
problems. Weighted spaces for the analysis of elliptic operators on non-smooth domains 
have e.g. been used in [28]. Weighted spaces have also used to analyze degenerate par- 
abolic equations, see e.g. [14, 27, 21]. Our analysis connects these two applications of 
weighted spaces. 

Let £ = (0, 1), £ = (0,°o) or £ = (-°o,0) and let d{x) = dist(x,(9£). For keN, our 
norms are given as the sum of two weighted Sobolev norms: 

IfUiE) - 4nf^_(N'+'5''+V+IL2(£) + ||^af+V-IL2(£)). (2.21) 

In particular, in the limit e — >^ the homogeneous norm in (2.21) turns from a norm of 
order 3/ -|- 1 to a norm of order 4/ -|- 1 . We furthermore set 

WfWhiE) - iif^liiEY (2.22) 

We recall Hardy's inequality which holds for all J3 7^ - i and all / e C~((0,oo)): 

II^^/IIl2((0H) ^ Cp\\^'^^'fM(o,-)y (2.23) 
In particular, for fixed £ > the second term on the right hand side of (2.21) is estimated 
by the first one. 

Our estimates require a generalization of the above norms to the case of fractional 
derivatives. This generalization will be done with help of the Mellin transform. This 
transform has been widely used for elliptic boundary problems on conical domains (e.g. 
[28]). For any / e C~((0, «>)), its Mellin transform / is 

-A ^/„^ dx _ f 5„ 



/(A) = / x-^f{x) — = / e-''"F{u) du. (2.24) 
Jo X Jm 

Here and in the following we will frequently use the variables u = Inx and F{u) = f{x). By 
(2.24), appUcation of the MelUn transform on / corresponds to application of the two-sided 

Laplace transform on F . It is easy to see that xfx{X) = A/(A) and f{X) = /(A -|- J3) 
for any j3 € M. Furthermore, Plancherel's identity holds 

ll^"''/llL2((0H,f) = lk~''"^llL^(M,d«) = II/IIl2(9{a=)3)- (2.25) 
The strip of convergence is the set of A G , j32) x M C C where the integrand in (2.24) 
is absolutely convergent. Note that if / is such that x^^^f and x^^f are in L^((0,°o), y) 
for some /3i < P2, then the strip of absolute convergence contains the interval (/3i,/32) as 
can be seen by applying Holder's inequality. For any /3 in the strip of convergence of /, 
the inverse Mellin transform of / is 

= L „^V(A)«S(A) = / e'"F{X)^{X), (2.26) 
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where the hne integral is taken in direction of increasing 3X. The definition (2.26) does 
not depend on the choice of j3 G (j3i , jSi) since / is analytic in the strip of convergence. 

We are ready to give a definition of the norm in terms of the Melhn transform. The 
definition of the norm by Mellin transform and the definition by (2.21) differ by a constants 
Q for integer k. In our notation we do not differentiate between the two definitions of 
the norms. This does not change the result since all out estimates depend on constants 
Ck- In order to apply the Melhn transform, we first need to subtract the boundary data. 
Consider / G C^{E), where £ = (0, 1) or £ = (0,°°), i.e. / vanishes for x — > oo. Let 
C G C"~([0,oo)) be a smooth cut-off function satisfying ^ = 1 in [0, |], C = in [3:°°) 
and IId'^i^II/,-. < Q for all integer k. Moreover, we assume that l^{x) + ^(1 —x) = \ for 
all X e [0,00). If £ = [0, 1] then we set /^(x) = /WCW and f{x) = /(I -x)C{x), in 
particular supp/^ supp/« C [0, 1] and /^(x) -x) = f{x){^{x) + ^{l-x))= /(x). 

We now define for A: > 

ll/llz|([0,l]) •= ll/^llx|((0,oo)) + ll/^llx|((0,oo))- '^^■^'^^ 

It remains to define the norm for E = {Q,°°): Given k >0, let n^. be the largest integer 
smaller than 3k — ^, i.e. ni^ = \3k — 5I • In particular, if A: G No then nk = 3k—l. Let Vf be 
the Taylor polynomial of order of / at x = (if % = — 1, then we choose Vf = 0). We 
decompose / = /i +/o, where /i = ^Vf and define 

ll/llx|((OH) = 11^/11 + II/oIIl^ + [/o]x|> (2-28) 

where || • || is any fixed polynomial norm, e.g. the ^^-norm of the coefficients. Here, the 
homogeneous norm [-J^i, A: > 0, is given by 

[/o]x| = m\"'^'^^'fo\\L^^x=3k-l) (2-29) 

with the notation jx ~ min{|A|, i}. The equivalence of these norms with the characteriza- 
tion (1.7) when k is an integer follows by application of (2.25) and by repeated application 
of Hardy's inequality; the proof is given in [26]. We define X| as completion of C°°{E) 
with respect to (2.28). The notation X^{E) indicates that the completion is taken in the 
subspace of C°°{E) where additionally / = on dE. Note that the trace of / is controlled 
inX^{E)\fk> g, in particular 7^ for > ^. Finally, we define as the completion 
of C~(£) with respect to (2.27)-(2.28). We will also use corresponding parabolic norms 
and spaces: Generically, the norms are defined for (?,x) G [0,oo) xE =:Q. For < e < 1, 
k G No, we define 

( Of 00, 

the corresponding spaces are called TX^, TX^, TXg, where as before the superscript '°' 
indicates that the function also vanishes at the boundary x = and the superscript 
denotes the space obtained by taking the closure of C~([0,°o) x (0,oo)). If the domain of 
integration is 2 = [0,oo)^ or £ = [0,°°), we sometimes omit the domain in the notation of 
space and norm, i.e. we write X| forXg ((0,oo)) etc. 

Note that the choice of norms (2.29) is supported by the investigation of the linear 
operator in [26]. In particular, by [26, Theorem 3.2] for A:, > 0, e G [0, 3(2^1 ^ J . we have 



(2.31) 
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Figure 2. Transformation from wedge to strip. 

2.3. Norms for the pressure. We introduce norms and spaces to control the pressure q 
(cf. (2.15)). Unlike the spaces X\, which can be expressed in both real variables and 
Mellin variables (cf (2.28)), our norms for the pressure can only be expressed in terms of 
Mellin variables. Roughly speaking, we apply the Melhn transform in the radial direction 
(with respect to the tip of the wedge), but not in the angular variable. The norm to control 
the pressure is L°° in the radial frequency variables and IT in angular variables. Using 
the L~-norm in angular variables enables us to obtain estimates which are optimal in e. 
A standard approach using an L^(L^) -norm in the pressure would not capture the optimal 
e-dependence which is needed for the convergence to the limit model. There are several 
technical difficulties connected with the fact that we have to take the supremum in v. One 
of them is that the norms cannot be expressed in terms of physical variables. Moreover, 
complex interpolation as is possible for the trace norm (see Lemma 5.2) cannot be directly 
used for the norms for the pressure. 

We first define the space on the wedge. We introduce a coordinate transform which 
maps the wedge onto an infinite strip: We define the new variables if^s) by 



x = e"cos(ev), . / 1^ _ /e"cos(ev) -ee"sin(ev) 

>'= ie«sin(ev), ''^ 1^ j " Ue"sin(ev) e"cos(ev) 



(2.32) 



For later reference, we note that dxdy = e-^"dudv and 

du = xd^ +ydy, = -eyd^ + ^xdy, (2.33) 

dx = e~"cos(ev)5„- je'"sm{ev)dy, ^dy = e~"sin(ev)^„ - ^e^"cos{ev)dv. 

The coordinate transform (2.32) can be understood as sequence of the two transformations 
{x,ey) = r(cos 0,sin0) and (r, 0) — (e", ev), see Fig. 2. Let jx — inf{i, Furthermore, 
for ^ € C~(/r\(0,0)), let. ^be the Laplace transform of q with respect to u, where q{u,v) — 
q{x,y). Suppose that q satisfies (2.15) with ^ = (the equation for the Unearized pressure). 
In the transformed variables, the equation for q has the form: 

X^q{X,v) + (^)2^vv(A,v) =0 in M X (0, 1) 

with boundary conditions ^(A, 1) = (A + 1)/(A + 1) and qv{X,0) = 0. This explicit ex- 
pression of q in Mellin transformed variables motivates the definition of our norms. For 

any multi-index a = (ai,a2) G K x No, we define |a| = ai -|-a2 and A" = X'^^{^dv)"^. 
For any e C^{K\{0,0)) and for £ > 0, we set 

[lo]yi = I II sup \e^'\'\^'-'^A«^i%\\y _s (2.34) 
\a\=3e ve(0,l) ^ i> y J 

For technical reasons that will be explained later, we will use these homogeneous norms in 
particular for £ G Suppose that / G X| and suppose that q is defined as in (2.15) with 

rj = fx. Since q = fxOndiK and since by the definition of X| the Taylor polynomial of / of 
order is well-defined, we expect to have control on the supremum norm for derivatives 
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of q up to order nk—l, where we recall that iik = [3A; — |] . Indeed, such an estimate is 
given in Lemma 5.5. 

For q e C~(/:), let Vq be the Taylor polynomial of q at (0,0) of order - 1 (if e N 
then nj. - 1 = 3fc - 2). Let C : be a cut-off function such that = C,{r) with C, = 1 

in [0, g], ^ = outside [0, \] and such that < ^ < L We decompose q := qo + qi with 
qi = ^Vq and define for A:, ^ e M, the norm 

M\yI = \\'P<i\\r+ sup teo]r|. (2.35) 

^<e<k 

Here, || • \\-p is any fixed polynomial norm, e.g. the £^-norm of the coefficients. We do not 
include the homogeneous norms with f < i since this would lead to 'negative derivatives' 

, oo, 

in the definition of the following norm, cf. (2.38). The spaces 7g and Yg are defined by 
completion with respect to functions C^{K) respectively C~(/r\(0,0)) as before. By the 
above considerations, the polynomial is uniquely defined and the norm is well-defined. 
The space Y^{Q.) and its norm are defined by localizing the above definitions (analogously 
as for the definition of Xg{[Q, 1])). The space TY^ and its norm are defined analogously to 
TX^, i.e. 

l<i-\- j<k 

Let us remark that we believe that all homogeneous norms in (6.13) for all realms [j,k] can 
be bounded by the two extremal homogeneous norms {£ = £ = k). However, the proof 
of this interpolation inequaUty does not seem to be straightforward, in particular since 
the analyticity of the expressions is destroyed by the supremum in v so that the theory of 
complex interpolation does not seem to apply directly. This is the reason, why we include 
the information about all the intermediate homogeneous norms into the definition (2.35). 
This is not necessary in the definition of the norm for the space Z| in (2.28) since there we 
have an interpolation result at hand (see Lemma 5.2). 

The space Z* describes the regularity of functions g = A^q for g e . In view of (2.34), 
this suggests to consider for real i€ [5, A:], the homogeneous norms of type 

Nzl = L II sup e^^l^l(i-)A«M'lfo|||,2(5R,=3£-2). (2-37) 

\a\=3e-2 vG(0,l) ^' 

where the sums are taken over a = (ai, ai) S K+ x No using the notation |a| = ai + ai. 
Note that [^0] z | < °° imphes 1 1 II lj' = o (x^^"'^/^) ) for all ^0 € Zg . Correspondingly, we 
say ^ G Zg if there is a polynomial Vg in x,y of order — 3 (if A: G N then nk — 2> ~2>k — A) 
such that ^0 := ^ — C.Vg G zf for some radial cut-off C, = ^(r) with ^ = 1 in [0, 1] and 
^ = in [2,00). The corresponding norm is given by 

ll^llzf = sup Mz|- (2.38) 

The spaces Zg and Zg are defined by completion as before. The corresponding space Zg(i2) 
for the droplet case and its norm are defined analogously as before by localizing the above 
definitions. Note that the minimal value ^ = | in (2.38) is chosen such that the exponent 
I a| = 3£ — 2 in definition (2.37) stays non-negative. 

2.4. Compatibility conditions. Higher regularity for our solution requires compatibility 
conditions on the initial data (for both Darcy flow and thin-film equation): Indeed, let 
/ e TXl be a solution of (2.16), respectively (2.19). Since / = at x = 0, it follows that 
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^t^fix=o = 0- This translates to a compatibility condition for the initial data. It is obtained 
by consecutively replacing the time derivatives in dff\x=o by the spatial operators A4 and 
using (2.16) resp. (2.19). The corresponding condition needs to be satisfied for the 
initial data: 

/in satisfies compatibility conditions ensuring ^//jx=o = at f = (2.39) 

for all < / < A:. For example, the condition foxk=\ corresponds to Ne (/in , /in ) = /in • 
An analogous compatibiUty condition needs to be satisfied for the linear evolution / + 
Ag/ = g. In this case, we need /in and g to satisfy for all < / < A: 

/in and g satisfy compatibility conditions ensuring djf^^^Q = at ? = 0. (2.40) 

3. Statement and discussion of the results 

3.1. Statement of results. We have the following main results: We have well-posedness 
for the Darcy flow with moving contact line. In the regime of lubrication approximation, 
we have convergence of the solutions towards solutions of the thin-film equation. Further- 
more as a consequence, we obtain well-posedness for the thin-fihn equation. 

Theorem 3.1 (Darcy flow). Let k>l be an integer, e e (0, ^(^S+i))- ^^PPose that /^ e 

Xg ' '^ satisfies (2.39) and suppose that ||/nll y*+'/2 < Oik for some (small) universal con- 

stant a^. Then there is a unique global in time solution G TX^ of (2.16)e with initial 
data /j^. Furthermore, 

ll/^llrz*+i + ll?^llrr*+i - ^kWfLW^Mn, (3.1) 

where p^ = p^ ^ Yg and : K ^ K^^ is the coordinate transform in Lemma 6.5; 
see also (2.28), (2.34) for the definition of the norms. The constant in (3.1) is universal, in 

particular it does not depend on £. 

The above well-posedness result can also be stated for the Darcy flow in terms of the 
original variables: Suppose that the assumptions of Theorem 3.1 hold. Then there exists a 
unique classical solution of (1.2). In particular, if /^ is sufficiently smaU then h^ defined 
by dxhf^ = 1 + /i^ satisfies > for x > 0. 

We also have convergence for solutions of the Darcy flow to solutions of the thin-film 
equation. Furthermore, as suggested by the asymptotic expansion, in the limit £ — )• 0, the 
pressure p is independent of the vertical direction: 

Theorem 3.2 (Convergence). Suppose that the assumptions of Theorem 3.1 are satisfied. 
Let f^ be the solution of (2.16)e with initial data /^ and let p^ be the corresponding 

pressure. Suppose that H/j^ — /in|L*+i/2 -^0 as e ^ Ofor some fm S X^'^^^^. Then there 

exist f, p and a subsequence Ej — )• such that 

Wf^' - fWrxl and \\p^^ - p\\j.yk -^0 as j^oo, (3.2) 

where p^=p^oW^ £ with : K ^ K^' defined in (6.25). Furthermore, f e rxg+i 
solves (2.19) with initial data /n. The limit pressure p does not depend on the vertical 
direction, i.e. p = p{t,x). 

As a consequence of Theorem 3.2, the velocity field U = (y,W) in the limit e = 
is horizontal and does not depend on y, i.e. U = {V{t,x),0). By Theorem 3.2 and by 
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Proposition 5.1(2), the solutions converge also in terms of Sobolev norms. We e.g. have 

II/'-/IIl2c|*~1(k2)^0' l|V/^'-V?llLrc?-2(R+xjf)^0 ase^O. 
For the case of a droplet as initial data, we have short-time existence: 

Theorem 3.3 (Droplet). Let k> I be an integer, E = (0, 1). Suppose that hf^ G H^{E) 
with hf^ >Oin{—l, 1) and with hf^ = 0, \hf^J = lon dE. Suppose that f^^ := [hf^ — 5(1 — 

x^)]x € Xg^^^^{E) satisfies the compatibility condition (2.39). Then there is a time T > 
such that for every £ G (0, 3(2^+1) )' ^^^^^ ^ unique short-time solution h^ of (1.2) with 
initial data hf^ (where h^ describes the profile of the propagating liquid). Furthermore, 
f := [h^ - 1(1 -x^)]^ G TXI+\{0,t) X E) satisfies 

ll/^llr4+l((0,T)x£) ^ Ck\\ftn\\xM/2^j,y (3.3) 

The solution depends continuously on the initial data. Furthermore "/llrx* a* 
j °°for a subsequence Sj — )• and f solves (2.19). 

Theorem 3.3 also shows that any solution immediately assumes a regularity of order ^ 
where we recall that e is related to the opening angle. This is the maximal regvilarity which 
can be expected in a non-smooth domain with opening angle of order e, see [28]. 

K+ ' 

Corollary 3.4. Any solution of f as in Theorem 3.3 satisfies / G Xg ^ (£) for any fixed 
positive time, where K is the largest integer such that K < 5 (jl — !)• 

Indeed, this follows by a bootstrap argument using (3.3): If fin G Xg'^^^'^{E), then we 
also have L?-{X^+^ ((0, t) x £)) which yields / G X|+i (E) for almost every fixed positive 

time to > 0. Application of Theorem 3.3 then implies / G L^(Xg^^^^((fo, t) x E)) for all 
t > to. Now, for any 5, we may repeatedly apply this argument for time steps of size ^ 
which yields the assertion of Corollary 3.4. 

Our analysis also yields the following new existence, uniqueness and regularity result 
for classical solutions of the thin-film equation: 

Theorem 3.5 (Thin-film equation). Let k>\ be integer 

(1) There is > such that for any fa G X^'^^^'^ with ||/in|| k+i/i < Ol.k and such that 

(2.39) is satisfied, there is a unique global in time solution f G TX^ of (2.19) with 
initial data fin- Furthermore, < Q||/in|Lj:+i/2. 

Aq 

(2) Let fin G Xq"*"^^^ and suppose that the analogous assumptions as in Theorem 3.3 
holds. Then there is a short time solution f G TX^^ of (2.19). Furthermore, 

ll/ellr4+i((o,T)x£) ^ ^k\\fAxl-^y^(Ey 

Note that the existence for weak solutions of the thin-film equation (1 .4) in the complete 
wetting regime (zero contact angle) is well understood, see e.g. [7, 8]; uniqueness of weak 
solutions is still an open problem. Existence and uniqueness of classical solutions has been 
shown in [21, 20]. All the above results address the case of complete wetting where the 
liquid attains a zero contact angle at the triple point. There are only few results for the 
partial wetting regime where existence (but not uniqueness) of weak solutions is proved 
[34, 9]. Well-posedness for classical solutions with non-zero contact angle for a related 
model has been shown in [25]. In this paper, we give the first existence and uniqueness 
result for (1 .4) in the partial wetting regime. We hope that the techniques developed in this 
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paper can also be applied to more complicated systems such as the Stokes flow with various 
boundary conditions at the liquid-solid interface or for fluid models where the contact angle 
condition at the triple point is different. 

In the following, we do not explicitly write A:-dependence of constants, i.e. we write 
C = Q. 

3.2. Formal lubrication approximation. We formally show how the Darcy flow (2.16) 
converges to the thin-film equation (2.19). For this, we show convergence of both linear 
and nonUnear operator in (2.16), i.e. Ag Aq and Ng ^ No as e ^ 0. The argument is 
based on an asymptotic expansion of the (e-dependent) pressure p'^ in ey (cf. (2.8)): 

p^t,x,y) =poit,x)+ E^P2{t,x, I) + 0{£'^). 

Our aim is to solve (2.8) up to first order in e, i.e. 

^.V + (H)V = ini^, ,,,, 

^ir = P,Vo = 0- 
Indeed, the solution of (3.4) has the asymptotic expansion 

p'{t,x,y) - /,_^(/-/,2)/„, + 0(e4). 
Inserting this asymptotic expansion into (2.8), we obtain 

BtUt,x) = B%{t,x) -hj^-hf^ + 0{e'') = -{hf^), + 0{e^). (3.5) 

where h = x + /q <j!>. The asymptotic expression of the Unear operator Ag follows as a 
special case of (3.5) by setting ^ = or equivalently h=x: 

Aef ^'=i'* -(fig/.). ^'='^ {xf^)^ + 0{e^) "-'^'^ Aof + O{e% 

which implies Ag — > Aq. The convergence Ne No can be seen similarly: With the notation 
3> = Jq (p dx = h — X, we have 



Ne{(p,f) ^='^ (pxf^\x=0 + {Btfx)x-{Befx)x + O{e^) 
^xxfxx\x=0 - {hfxx)xx + {Xfxx)xx + O(e') 



p=0 ~ {hfxx)xx + {xfxx). 

which formally proves the convergence Ne^ Nq. 

4. Proof of Theorems 3.1-3.5 

We give the proof of the Theorems 3. 1-3.5 already in this section so that the reader may 
get an overview of the structure of the proof. Some parts of the proof are based on results 
and estimates which are given in detail in the later part of this work. 

4.1. Proof of Theorem 3.1. The proof of Theorem 3.1 proceeds by an application of a 
contraction principle. It is based on maximal regularity for the linear operator Ag and 
corresponding bounds for the operator A^g. In the following we drop the superscript e in 
the notation if the meaning is clear from the context, e.g. f = and correspondingly for 
the other functions used. The maximal regularity estimate is the foUowing: 
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Proposition 4.1. Let k>\ be an integer, e e (0, jpfpyj)- Suppose that G and 
g e TXg~^ satisfy (2.40). Then there is a unique global in time solution e TX^ of 

Ut+Aef = g /orxe(0,-), 

\/ = forx = 0. ^ ■ ^ 

with initial data /;„. Furthermore, for some uniform constant C > and allx>Q and with 

Qx — (0,t) X (0,oo), we have 

ll/llrxi+Hex) - ^ (ll^llrx|(er) + ll^lli2(x«(Gr)) + ll/i"llxf'/2((o.oo)) + il/llc!'L2(e,))) (4.2) 

The estimate on the nonlinear operator is stated in the following proposition: 

Proposition 4.2. Suppose that the assumptions of Theorem 3. 1 are satisfied. Suppose for 
i = 1,2 thatf e with WfiW^^k+i. Then 

l|A^a/l,/l)-A^e(/2,/2)|lrx|(e.) + ll^^(/i'/l)-^^(/2'/2)llL2(xO(e.)) 

< C'll/i -/2|lrz|+i(eT) (11-^1 llrx|+'(et) + \\f^\\Txl+\Qr))- ("^-^^ 
Proposition 4.1 has been derived in [26J. Note the extra term II/IIc<'l2(2t)) o° the right- 
hand side of (4.2) which is due to the slightly different definition of the norms ||/|| j^^t+i (g^) • 
This extra term is needed for the nonlinear estimate (4.3). The proof of Proposition 4.2 is 
shown in Section 6. The estimates in the above two propositions hold for every chosen 
time interval. The constants do not depend on this interval. 

The proof of Theorem 3.1 now follows by application of a contraction argument: For 
5 > to be fixed later, we set 

£ := {/er4 : ll/ll^^l < 5}. (4.4) 

The operator Ss^f) is defined for any f E E as the solution of (2.16) with fixed initial 
data /in and right hand side Ne{f,f). Let /i,/2 G E, with the same initial data and let 
/ := fi — f2- In particular, 5e(/i) —S(;{f2) solves (2.16) with vanishing initial data and 
right hand side Ng{fi,fi) — Ng {fi , /2 ) • By standard interpolation, we have for any f E E, 

11/11 

Hence, there is a small but universal constant T such that for t = To, we can absorb the 
last term on the right hand side of (4.2) (increasing the constant in the estimate by some 
universal factor) to get 

ll-^llrx|+i{e.) ^ (ll^ll(rx|nL2{xO)(ex) + ll/inll4+'/2((o,oo)))' ^^.6) 

where we use the notation || • || jXjnL2(z|') ■~ II ' ll^^c + II ' \\l^{x^)- ^'^■^^ ^ ^^^^ 
(4.3), we hence get 

\\Se{fl)-Se{f2)\\rxt^ < C\\Ne{fufy)-Ne{f2,f2)\\ 
(4.3) 

< C'll/i -/2|lrx*+' (ll/illrxl+i + ll/2|lrx-|+0 

< C5||/i-/2||y^.+,. 

Hence, Sg is a contraction if 5 > and T are chosen sufficiently small. Similarly, by (4.6), 
(4.3) and since A^e(O) = 0, we get 

\Mf)\\rxt' C||iV,(/)||^^,nL2(xO)+C||/in||^.+i/2 Ca+C5^ 
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and hence S£{E) C E for 5 and a = (x{5) sufficiently small. Therefore, application of 
Banach's Fix-point Theorem yields existence and uniqueness of a solution of (2.16) on the 
time interval (0, 1). In order to recover long-time existence, we use dissipation of energy 
(1.5). Indeed, by (1.5), we have 

\\f\hLUQr) ^ C||/in|L2((o,.)) < C||/i„||^,^i/2. (4.7) 

Using this estimate instead of (4.5), we get estimates that are independent of the time 
interval. This shows long-time existence and also the estimate of / in (3.1). In order to 
conclude the proof of Theorem 3.1, it remains to prove the uniform bound (3.1) on p. This 
estimate follows from the estimate in Proposition 6.7. 

4.2. Proof of Theorem 3.2. By Theorem 3.1, we have the uniform bound ll/^^H^j^t+i < 

Ca for all e > 0. We use the optimal decomposition into high and low frequencies =^ 
f+ +/- from (5.6). By (5.6), this decomposition commutes with the time derivative. We 

have uniform bounds on the norms ||-^^-'^'^/^;^^''^V+IIl2((o,oo)2) ™d ||-^^/^i''^V-llL2((0H^) 
for all i, j with i + j <k+l. Standard compactness appUed to both /| then show that in 
particular there is f = f^+f^ G TX^ and a subsequence Ej — )• such that 

Wf^-fWrxl^O as7^~. (4.8) 

Now, let p^', p^i be the pressure related to and f^i. By (6.35), we then have 

\\f'-f'\\TYi < CWr-Mrxl {\\f%xl +ll/'^llrz|) CaWr-f'Wjxl ^0 

; ; 7 i 

as i, 7 — > °°, 7 < i. This shows that p^> converges in TY^ thus concluding the proof of (3.2). 

It remains to show that / solves the thin-film equation (2.19). By the above uniform 
bounds, by (5.2) and by (5.11), we have ^ / in L^{H^), ff ^ /, in and p^ ^ p 

va.l}{H^) with ||p^||£2(£2) < Ce-J> 0. The boundary condition p^ =hxx{i + £'^hly^ hence 
implies that in the limit e = 0, we get p = hxx- We will show that in the Umit e = 0, is a 
solution of the thin-film equation, where we recall that ~h^ — l. In particular, ^ hin 
L^{H'^) and ht in I?{L?-). Correspondingly, we also have convergence of the velocity 

= (v^,w^) = Vep^ u = (v,w) = {px,0) in L?{L^). The transition to the Umit now 
can be conveniently done in terms of the continuity equation: By conservation of mass for 
(2.16), we have 

hf+[J v^dxj = 0. (4.9) 
In the limit £ ^ and in view of the above discussion, (4.9) turns into 

h, + ( vdx^^ = ht + {hpx)x = h, + {hhxxx)x = 0. (4.10) 

4.3 . Proof of Theorem 3.3. We prove this theorem by an application of the Inverse Func- 
tion Theorem. For this, we linearize Ce at an 'approximate solution' w, constructed with 
the help of an extension lemma. We then show boundedness and differentiability for £e 
and invertibility and maximal regularity for its linearization 5£g at w. We keep the details 
brief and refer to similar arguments in [14, 20, 3]. 

We define the 'spatial part' Ae of the operator by £ef = Dft+Asf, i.e. 

Aef := D,f-dx\{{l-x)s-+xs+){f+{x{l-x))x)+B{(- 

L ^{l+£^if +{x(l-x))x)^)^ 

see (2.9). We will use the following extension Lenoma 
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Lemma 4.3 (Extension Lemma). Let k^fi with k>\. For any /in € x'^'^^^^{E), gi^ G 
Xg ^^^{E) there exists w e 7'Z|+^([0,<») x E) such thatw^^^o = fm, = gj^ and 

Il^llr4+l([0.~)x£) - C"(ll/mllx*+l/2(£) + |km||j,*-l/2(^^). (4.11) 

The extension can be constructed by gluing together solutions of the linear equation 
given by Proposition 4. 1 . The methods used in [20] can also be used for our equation so 
that we will not present the argument here. For /n € X|, let gin = — ^e/in and we choose 
w G rx|+^([0,oo) X E) as in Lemma 4.3. In particular, 

^ew\t=Q = Ddtw\,=Q + A£f in ^''= ' gin-^in = (4.12) 

and w may in this sense be called an approximate solution. Let 5£g be the linearization 
of £e around w. We have boundedness and differentiabiUty of Ce and boundedness of 
{8M(w))~^ for T small enough: 

Proposition 4.4. Let k> 1, k eN and suppose that fm G x'^^^^^{E) with f fmdx = and 
g e TX^{Q^) satisfy (2.40). Let w € X|+i([0,°o) x E) be defined as in Lemma 4.3 with 
gin := ^|r=o- Then for sufficiently small T > and with the notation Qx = (0, t) x E, there 
exists a unique f e rx|+' {Qx), solution of5jCe{w)f = g in Qx with f = 0on (0, t) x dE 
and with initial data f = fin- Furthermore, 

ll/llrx|+i(&) ^ ^(ll^llrx|(&) + ll/in|l4+i/2(£)). (4.13) 

Note that the condition f{x)dx = is preserved by the flow generated by 5Ce{w)f . 
We also have differentiability of Ce in a neighborhood of w: 

Proposition 4.5. Suppose that fin € X^'^'^^^{E) satisfies (2.39) and J^findx = and let w 
be defined as in Lemma 4.3 with gin = —Aefn- Then for sufficiently small T > there is 
a > such that M.g, : TX^'^^{Qx) Xg'^^^^{E) x TXg{Qx) is bounded and continuously 
differentiable in the a-neighborhood ofw in TXg'^^{Qx). 

The proof of the above two propositions is given in Section 7. Using the above two 
propositions, the proof of Theorem 3.3 follows by application of the inverse function the- 
orem: We claim that the operator 

Me:TX^+\Q,)^TXl+^'^{E)xTXl{Q,) with / ^ Me{f) := {f\t=oXef) 

is bounded, continuously differentiable near w and 5A^£(w) is invertible with bounded 
inverse for T small enough. We define v := Ceiw). By the inverse mapping theorem there 
is a neighborhood of w and a neighborhood of (/n, v) where Me is a diffeomorphism. By 
(4.12) we have v\t=o = 0. Since C^{Qx) is dense in TX^{Qx), it foUows that IkllrxKe^) ^ 
for T ^ 0. Hence, there is f <5 (0, t) and a function v GY^ with v = for f G (0, f) and 
such that {0,v) is sufficiently near (0, v). Hence, there is / G with Me{f ) = (0, v). The 
function / is a solution of Cgf = and hence h{x) = x{l —x) + fQf{x')dx' is a solution of 
(2.1 1) for t e (0, f ), thus concluding the proof of Theorem 3.3. 

4.4. Proof of Theorem 3.5. By Theorem 3.2, we have existence and regularity of so- 
lutions of (2.11) for initial data which are close to the infinite wedge. In order to show 
Theorem 3.5(1) it hence remains to prove uniqueness of solutions. For this, it is enough to 
show that the corresponding results in Propositions 4.1 and 4.2 also hold in the case e = 
and for the operators Aq and A^o- The estimate for existence, regularity and uniqueness in 
the case of half-space then foUows by the same fix-point argument as for Theorem 3.1. 
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The e = version of Proposition 4.1 has been proved in [26]. The £ — 0-version of 
Proposition 4.2 can be obtained by analogous estimates as the one's applied in the proof of 
Proposition 8.1 in [21]. In fact, the estimate is easier in our situation since we only need 
for an estimate in weighted Sobolev spaces, not the interpolation spaces used in [21]. Fi- 
nally, we note that the local result in Theorem 3.5 can be obtained by standard localization 
techniques. The argument can be performed analogously as our locahzation argument in 
Section 7; the argument is easier since for e = 0, the norms || • \\^k are local. We also refer 
to a similar localization argument in [20], performed for a thin-film equation in a Holder 
space setting. 

5. Estimates in weighted spaces 
We recall some basic properties of the space X| (see Proposition 2.3 of [26]): 

Proposition 5.1. Let A: e M with k>0 and let f e X\. 

(1) For < e < e', we have X| C X* and for all f e X| 

\f\x\ < (5-1) 

(2) For any £i,£iG No with 0<£i<3k and < £2 < 3A: - j, we have 

\\d!'f\\L2 + \\dtf\\L- < C\\f\\^,. (5.2) 

(3) /// G 4, then dj e xl-''\ d^f, d^f G and 

ll/.||^*-i/2 + WfxAlxr + ^ (5.3) 

The norm || • controls a scale of weighted Sobolev spaces of fractional order: 
Lemma 5.2. Let e G [0, 1 ), A > and fo G X^- Then for every <£ <k,we have 

||A^^+VVo!lz.2(CR;i=3^_l) + ||/o|li2(9{;L=3^_l) <C||/o||;f|. (5.4) 

We have the following characterization of the homogeneous norms: 
Lemma 5.3. LetkG Nq, e G [0, 3(2^+1) )■ Then for fo GXg we have 

c\fo]xl< ^Jnf^^ {\\^''^'U\\LH^^^ <C[/o]4, (5.5) 

where Pk = 3k — ^. Up to multiplication by a constant that only depends on k, for any even 
integer M with M>kthe minimum in the right hand side of (5.5) is achieved by 

= (l-(itan)^(eA))/o(A), and /L(A) = (itan)^(eA)/o(A). (5.6) 

Furthermore, 

\X\'\U{X)\ < Cm*|/(A)|, {\f\%{X)\ < Cm'|/(A)|. (5.7) 

The proof of Proposition 5.1, Lemma 5.3 and Lemma 5.2 is given in [26]. 

The space X| is an algebra as is proved below. Note that there is a proof for the algebra 
property in [26]). The advantage of the result below is that it also applies to the case k=\, 
also the proof is simpler than in [26]: 

Lemma 5.4. For A: G N with k>\ and f,gG X|, we have fg G X| and 

WfgWxl < cWfhiMxt (5-8) 
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Proof. We decompose / = /i +/o and g = gi+go where f\ —VfC^gi =^gC and where 
'Pf,Vg are the Taylor polynomials of /, g of order ?>k—\. The cut-off function C, is defined 
as for the definition of (2.28). In particular, /o,^o G We need to estimate the products 
/i5i,/o^b/i5oand/ogo- Clearly ||/i^i||;f| < C||/i||j^|||gi||^|. As in (5.6), we decompose 
/o = /+ + /- and ^0 = ^+ + ^- • The mixed term figo is estimated as follows 

[/igolxl < C||x^+^54^(/ig+)||^2((0H) + ll(D'^^''(/i^-)llL2({0H) 

< c||/||^|(||x*+'54VIIl2 + ||x5VIIz.O < 

The second estimate follows since all derivatives of f\ are supported in (g , |) and since x 
is of order 1 in this interval. The estimate of fygi proceeds analogously. 

It remains to show the estimate for the product /o^o = f+g+ + .f+g- + .f-g+ +f-g-- In 
particular, since fo,go &X^, we have |/o|, |^o| ^ Cx^*"' for small x. It follows that |/o^o| ^ 
^^6*-2 ^ (7x^^-2 (jjjjQg ^ > ^-j^ jjj particular fygo e We show the estimate for the term 
high-high-frequency product f-g-- Recall that the Mellin transform for the product of 
functions can be expressed as a convolution for the Melhn transformed functions. Hence, 

\f-8-\xl ^='' P''+V'/^^^/-(A-T,)|_(T,)^(T,)||^,(3,;^^3,_I), (5.9) 

where n = minj^, and where a e M is chosen such that the product /-(A — ■)§-{■) 
is absolutely integrable on the line 9^7] = a. Since ^ < |, the right-hand side of (5.9) can 
be estimated by replacing /x*^ by (^)*. Using the binomial formula A^*"*"^ = X),c;fe,Tj'(A — 
j^^3i;+i-!^ (5.9) is bounded by above by the sum of the terms 



hAf-,g-)=l , / (i)^(A-T,)3^+l-'/_(A-T,)T,'|r(T,)^S(T,) 



where < i < 3A: + 1 . By symmetry, it is enough to estimate the terms with i< + j . Note 
that the integrand of the inner integral above is analytic as a function of t/ . In particular, the 
value of the integral does not depend on a. This argument works since we have avoided to 
replace A by | A | in our proof. By Young's inequality for convolutions and by the Cauchy- 
Schwarz inequahty, we have 

\\F*G\\^H%x=fi) < Qll(l + l^l^)^llL2(9?;i=/3,)l|G|L2(9;A=fc): (5-10) 
which holds for all 5 > j. With the choice a = i and 5 = 1, we get 

C|!M'-^|Ap'=-'-+7o||,2(9,,.3.-,--i)liM^(l + |Ar'+')fo|L2(95,.,^ 

The last estimate follows from (5.4) using < / < jk+ j (and since k > 1). The estimate 
of the terms f+g-, f-g+ and f+g+ proceeds analogously (see also the related proof in 
[26]). This concludes the proof of the Lemma. □ 

Also the space is embedded into classical Sobolev spaces: For any multi-index a e 
]R>o X No, we set |a| = ai -1- 052. For a e we also usethenotationZ)^ = 5"'(i^j,)"2. 
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Lemma 5.5. Letk^M. with A > 1, e e (0, 1) and let q e Fg. Then for all £i, li G No with 
3 

2 



l\ < 3A: — 1 and £2 < 3k — I (or equivalently ^2 < «fc — IJ. have 



I l|OsMlzrLi(jf)+ I sup ||DgMU.(j,) < Clklly*, (5.11) 

0<|a|<£i 0<|a|<^2y^('''^) 

where the sums are taken over multiindices a G I?. 

00, 

Proof of Lemma 5.5. Let q = qo + q\ with <7o <= i'g be the decomposition as in (2.34). By 
definition, qi solves (5.11). It suffices to show the estimate (5.11) for ^0 S Ye- apphca- 
tion of Plancherel's identity, we have 

L sup ||£'£^o||l2((o,oo)) ^ ^ L ll'5«"'(i'5v)"'9ollz.2((oH) 
|a|=3£-iye(0,A:) |a|=3£-ive(0,l) 

(2.34) 

< C 21 sup ||AVo|li2(m=3£-3) < Cfeolvt- 

|ce|=3^-lvG(0,l) ^' 

This yields the I? estimate, the supremum estimate follows by standard interpolation. □ 

For any / e its trace at ? = is well-defined in xl'^^'^: 

Lemma 5.6. Let k,Y& Nq. For f e TX^+^ we have 

\\d'f\\co^^^.m^ < c||a/+VlL2(^|)+c||a//||^,(^..:). (5.12) 

Proof. It suffices to give the proof for ; — (for ; > consider F = dlf instead). By 
an approximation argument it is enough to consider / S C~([0,oo)^) . We decompose 
/ = /o +/i where /i —VfC, and where Vf is the Taylor polynomial of order n<:+i/2 = 3k of 
/ at X = 0; in particular, /o € ' . In order to avoid fractional derivatives which appear 
in the definition of the norm forXg*+i/2 , we use the equivalence c[/o] *+i/2 < {Afo,fo)x'' ^ 
C[fo] k+i/2 which holds for all /o e Xg'^^^^. A proof of this equivalence is given for A: = 
in [26, Lemma 4.6], the argument used there also applies for general e N. The estimate 
of the homogeneous part is then easy: For f,g e Xg^'^^ nx|+\ we have (Ag/o,^o)x* — 
(/o,Ag^o)x*- Hence, integrating in time from infinity (where / = 0), we obtain 



sup[/o]t*+i/2 < Csup(Ag/o,/o)x* < C I (Ag/o,/of)x* 

t t ^ Jo " 



dt 



< 



C[Aefo\L2{xl)lfot]L2{Xl) ^ <^[/0]l2(x|+1)[/0?]z,2(x|)- 



It remains to give the corresponding estimate for f\ : That is, we need to estimate the 
coefficients of the Taylor polynomial Vf. We show the estimate for the highest order 
coefficient of Vf. Up to a constant it is given by F{0) where F := d^'^f. We extend F 
symmetrically as an even function defined for all e K by setting F{t, —x) := F{t,x). 
We claim that 

sup|F(0)| < C(||F,||^2(K2) + ||F;^|L2(K2) + ||F||^2(K2))- (5.13) 
Indeed, (5.13) can be derived by taking the Fourier transform F(t7 , ^ ) of F(?,x) 

sup|f(«)| < JJlFi < (//(rH^,|^)"'(//(> + lnl^+l«ftl^p)"'. 



20 



HANS KNUPFER AND NADER MASMOUDI 



The first integral on the right hand side is bounded: One the one hand we have 

1 



dt,dr] < C < oo. 



(5.14) 

(5.15) 

□ 



/|||<lor|7,|<l 1 + |t?|2 + |^|6 

Also, with the coordinate transform §^ = tj^A^ and =r\^l^dX 

r r I 17 ^ 1,1. d^dr^ < C Ttj-^/^ r \ dUr] <C < 
Ji Ji \ri\^ + \^\^ ^ ' - Ji ' Jo 1 + |A|6 ' - 

This concludes the proof of (5.13) and thus of the lenrnia. 

6. Uniform estimates for the operator in the half-space 
6.1. Linear pressure estimates. We derive estimates for the pressure p <eY^ 

Proposition 6.1. Let e M with > 0, e € (0, 3(2^+1) )• Then for any / G X| and any 

g there is a unique solution p e of 



Furthermore, we have 



in K, 
on diK, 
on doK. 



llPllrf < + 11/11.*) • 




(6.1) 



(6.2) 



We first address the situation of homogeneous data f GXg and g G Z|: The coordinate 
transform (2.32) leads to the following model for p understood as a fimction of (m, v): 

d^P + {-ed,)^P = e^"g in MX (0,1), 

p{u,-) = e-^fuiu) forv=l, (6.3) 

Pv{u,-) =0 forv = 0. 

AppUcation of the Laplace transform (2.24) in terms of u yields 



X^p{X,v) + ad,)^p{X,v) =g{X-2,v) inKx(0,l), 



/5(A,.) = (A + 1)/(A + 1) 

P,{X,-) 

Explicit solution of (6.4) for g = yields: 



for V = 1, 
for V = 0. 



(6.4) 



Lemma 6.2. Suppose that the assumptions of Proposition 6.1 hold. Suppose that / € X* 
and g = 0. Then there is a unique solution q '■= p & Yg of (6.1). It is given by 



Moreover, \\q\\yk < C\\f\\^k. 

Proof Clearly (6.5) solves (6.4). The estimate WpWyk < C\\f\\fk follows from 



(6.5) 



sin(eAv) 



cos(eA) 



< 



cos(eAv) 



cos(eA) 



(6.6) 



< Ce^l^K"-!) < Ce'^^m^-^) 

for V e (0,1) which also implies uniqueness. Indeed, for |eA| < 1, we have n = |A|, 

aid 

□ 



^^^^1 < C and |ele^l(^-i)| > e"!. For |eA| > 1, we have |cos(eAv)| < 4e\'^\' and 



I cos(eA) 

I sin(eA)| > igl^^l. This proves (6.6). See also [26]. 



Explicitly solving (6.3) for / = yields: 
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Lemma 6.3. Suppose that the assumptions of Proposition 6.1 hold. Suppose that / = 
and g e Zg. Then there is a unique solution w:= p GYg of (6.1). It is given by 

e 



+ 



Acos(eA) 



sin(eA(v-l)) / cos{eXz)g{X -2,z)dz. 
Jo 



(6.7) 



Moreover, WwWyk < CWgW^k 



Proof. The solution w can be expressed in terms of the Green function G{v,z) by 

r e^X^Giv,z) + d^G{v,z) = 5,=, inH, 
w(v) = / G(v,z)e2|(A-2,z)Jz, where ^ G = forv=l, 

[ Gv = forv = 0. 

Since G is harmonic away from v = z and continuous at v = z, it must be of the form 



G(v,z) = 



Csin(eA(z- l))cos(eAv) if v < z, 
Ccos(eAz)sin(eA(v— 1)) if z < v, 



for some constant C to be determined. Taking the derivative of the above equation in v, we 
get 

-J ^ _ / -CeAsin(eA(z- l))sin(eAs) if v<z, 
c/yU(v,z) - j ceAcos(eA(v-l))cos(eAz) if z<v. 

Since the jimip of dyG{v,z) at v = z is 1, namely [5vG(v,z)] = 1, we deduce that 
1 



C 



cos(eA(z— l))cos(eAz) + sin(eA(z- l))sin(eAz) = eAcos(eA) 



which impUes (6.7). We next give the proof of the estimate: We will use that 



|sin(eA(v-l))| <Cel^^l(i-''), |cos(eAv)| < Cgl^^l" and 
Hence, using the notation g{a) = g{X — 2, a), we infer that 



cos(eAv) 



cos(eA) 



<Cele^l(''-i). 



[\m{eX{a-l))g{a)da < C /^I^^Ki-^^e-T^^i"^' sup \e^-^^^-''^g{a)\da 



iv 

C \eM 

eA 



(l-v) 



sup \e 

(T€(0,1) 



(l-CT) 



ml 



(6.8) 



and 



f\os{eXa)g{a)da < C f e\^M<ye-^^{i-<y) ^up \e^^('-''^g{a)\da 

Jo ■^O CTG(0,1) 



- eA 



I Ml 
sup \e 2 



(i-a) 



(6.9) 



CT6(0,1) 

By (6.8)-(6.9) and in view of (6.7), we deduce that 



M>(A,v)| < sup |e^(i-^)^(A-2,v)|. 



|A| 



(6.10) 



v€(0,l) 
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We calculate the first derivative, 

,?2 



dyw{^,v) = sin(£Av) / sin(eA(v— l))g(A — 2,v)i/v 

g 

- -^cos(eAv)sin(eA(v- 1))|(A -2,v) 

^ COS( £^ ) 
a2 



+ 



-cos(eA(v— 1)) / cos(eAv)g(A — 2,v)dv 



cos(eA) 

^ sin(eA(v-l))cos(eAv)|(A-2,v). (6.11) 



Acos(eA) 

A similar calculation as before shows that 

|^Wv(A,v)| < ^e^^'-'^ sup |e^(i-^)^(A-2,v)|. (6.12) 

v€(0,l) 

Multiplication of both sides of (6.10)-(6.12) by |A|*^ yields higher regularity in the radial 
variables. Higher regularity in v follows by (6.10), (6.12) and repeated application of 

|A^(i)2wvv(A,v)| |A^+2m)(A,v)| + |A*^(A-2,v)|. (6.13) 
Estimates (6.10), (6.12) and (6.13) imply 

^ sup |eT^(i-")A3*+Vw(A,v)| < C ^ sup \e^^^-''^A^''n''g{?i -2,v)\. 

|a|=3it+2vG(0,l) |a|=3it''G(0,l) 

Estimate (6.2) follows by taking the L^-norm on the line 3iX = 3fe — j on both sides. □ 

Proof of Proposition 6.1. Let Vf be the Taylor polynomial of / at x = of order = 
3 A: — 1; let Vg be the Taylor polynomial of g at {x,y) = of order ni^ — S. Let Vp be the 
polynomial solving (6.1) with g and / replaced by Vg and Vf. Existence and uniqueness 
of this polynomial solution follows from a straightforward calculation and furthermore 

WVpWr < C{\\Vf\\r + \\Vg\\r), (6.14) 

see also the proof of Lemma 4.2 in [26] as well as the proof of Lemma 6.4 below. Let 

C : -> K be a cut-off such that C = C('') with C = 1 in [0, g], C = outside [0, i] and 
such that < C < 1 and let (0, M be given by ^{x) = l{\{x,x)\). We define 

Pi ■■= 'PpC, gi ■= ^ePi, /i := Vpdx, (6.15) 

00 T oo, 

and ^0 ^ — ^1 € Zg and /o := f — fi G Xg. Furthermore, let pQ be the solution of 
Apo = go in K with po = rjo on diK and dyPo = on ^o^- By the previous two lenomas, 
there exists such a solution satisfying 

IIpoIIk,* < C{\\gor£l + \\fo\\xl)- (6.16) 
Hence, p'.= po+ Pi is a solution of (6.1) and satisfies the desired estimate. □ 

Since e has the scaling of vertical length, one could expect that there is only uniform 
control on the norm j{py)\r- But it turns out that even (|)^(py)|r is bounded uniformly 
for e > 0. The proof of this statement is given in the next lemma: 

Lemma 6.4. Let k>l,s€{0, ^(^ik+i) )• ^^^'^ solution p € of (6.1) satisfies 

mPx)\r\\xl + i-efmPy)\r\\xl < C(||/|Ui + IklUi). (6.17) 



WELL-POSEDNESS OF DARCY'S FLOW AND LUBRICATION APPROXIMATION 



23 



Proof. Analogously as in the proof of Proposition 6.1, we decompose p = p\ + po where 
pi encodes the expansion at the boundary and where po € Ye^^ is the solution of (6.1) with 
corresponding homogeneous data /o € and G ■ Furthermore, let pQ = qQ + wq 
where is the solution of (6.1) with boundary data /o and with right hand side ^ = 0. 
Correspondingly, wq is the solution with right hand side and with boundary data / = 0. 
In the sequel, we present give the corresponding estimates to (6.17) for qo, wq and p\\ 
together these estimates imply (6.17). 

Estimate for qo: With the transformation (2.33), we need to show 

lk~^"^oHrllz| + (i)^lk~^"^o,vM|rllx| ^ C'||/o||^t+i. (6.18) 



Here, and in the following, by a slight abuse of notation we understand qo as a function of 



(m,v). By (6.5), since |sin(eA)/cos(eA)| < C;U and for |9^A| € (i 4A;- i) we obtain 



|e-2>,v„(A,l)| = \ek^ co"s(£(l 1 2)) + ^^-^"^^ + I - <^|£'^Wo(A+3)|. 

Multiplying this identity by |A and taking the L^-norm on the Une 9tA = 3A: - 5, 

we obtain the estimate for the second term on the left hand side of (6.18). The estimate for 
the first term proceeds analogously. 

Estimate for wq.- With the transform we need to show 

lk"'"v^0Hrllx| + (|)'lk"'">^0Hrllx| < C|ko|lzW- (6-19) 
Evaluating (6.1 1) at v = 1 we note that only the third term does not vanish, i.e. 

dMX,l) ^*'=' — / cos(eAv)|(A-2,v)dv. (6.20) 
cos(eA) Jo 

With the notation (p = (^)^e~^"wo,v«|r and since ^(A) = \{X+ 2)wo}{'k + 2, 1), we get 

cos('(A+2)) r-"-(^(^+^>")^»(^'-~)-^'^- ^'-^'^ 
For |eA| < 1, we have | < 1 and | cos(eA)| > i, |A| = and therefore 

|^(A)| < C|A| sup \go{X,v)\ < Cm sup |e^(i-^)|o(A,v)| (6.22) 

vG(o,i) ve(o,i) 

For \eX \ > 1, application of (6.9) yields 

|9(A)| < f sup |e^(i-^)^(A,v)| < Cm sup |e^(i-^)^(A,v)|, (6.23) 
v€(0,l) ve(o,i) 

where we also used that n = jin this case. The above two inequaUties together imply 

|A3^+V'9(A)| < C sup |e^(i-^)A3^+V'+'^o(A,v)| 
ve(o,i) 

for aU ^ > 0. Integrating the square of the above estimate on the line 9tA =3k—j yields 
[<PW = P''^V^^II,.(9t,=3.-.) < C|| sup k^(^-'')A3^+V^+V'o|||,2(9,,=3£_.) 

^ ^' v€(0,l) ^ ^' 



< C\\go\\ 



which concludes the estimate for the second term on the left hand side of (6.19). The 
estimate of the first term on the left hand side of (6.19) proceeds similarly. 
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Figure 3 . Coordinate transform onto moving domain 

Estimate for pi: Let Pp = T.i+j<eciijx'y^ be the polynomial which solves (6.1) with 
boundary data Vf and right hand side Vg, where Vf = Y^i^ix' is the Taylor polynomial of 
/ of order 3k — I and where Vg — Y,ij8ij^'y^ is the Taylor polynomial of g of order 3k — 4. 
Analogously as in the proof of Lemma 6.3, we need to show that the coefficients of the 
polynomial i^)^'Pp,y\r ^re bounded by the coefficients of /: Indeed, since by the condition 
we have d^Vp y^y^o = 0, it follows that a, i — for all / > 0. With the equation, i.e. 

d^d^,,^ = £2 {diVg^y - d^^Vp.y) , (6.24) 

we first get |a,;3| < e^HPgH and then iteratively |a;,27+i| < £^||'?'g|| for all i,j > 0. Fur- 
thermore since Vp^r = fx and again using (6.24), one can easily deduce that \ai^2j\ < 
Ce^{\bi+j\ + \\Vg\\) for all i,j > (we e.g. have 002 = e^{a20 + goo) and ao2 + ao2 = h). 
In particular, |ay | < e^dlP/Hp + ||'Pj||p) for all i>0,j> 1 which yields the desired esti- 
mate. □ 

6.2. Pull-back onto wedge. We need to measure the difference pi — p2, where pi is the 
solution for the pressure on the domain K^^ and p2 is the corresponding solution on K^^, 
see Proposition 6.7. For this, for given profile function (p, we introduce a pull-back from 
the perturbed wedge K''' to the unperturbed wedge K, see Fig. 3a). The estimates are 
nonlinear due to the geometry of the domain. 

Lemma 6.5. Let k €'N with k> I, e € {0, ^pfpyy)- Then for any (p e X| with W^W^k 
sufficiently small, there is a diffeomorphism W : K ^ Kl' of the form 

^'ix,y) = {x,y) + {OMx,y)) = {x,y). (6.25) 

Furthermore, analogously as in the definition (2.34), there is a decomposition Y=YO'^Vl 
such that \ffi = and where V is a polynomial of order 3k such that 

sup £ II sup |e^^l^l(i-)/+iA«£ljro|||^, 3,^1 < C\M^, (6.26) 

0<i<k\a\=ie+2 0<v<l ^' 

where we recall that n = inf { | A | , g }. Furthermore, we have 

\\eV.\r\\4 + \\Vy\r\\x^ < CM^k. (6.27) 

Proof. We construct Xj/ to be the solution of 

AeV^ = in^, 

\ir = Jofdx, on diK, (6.28) 
V/ = 0, on ^o^- 

As in the previous section, the argument is based on a decomposition the right hand side 
into a polynomial part and a homogeneous remainder. Since this decomposition proceeds 
analogously as in the proof of Proposition 6. 1 we only consider the case of homogeneous 
data. That is we assume that (p e Xg. With the transformation (2.32), (6.28) takes the form 
^uV+ i^d^fW = for (m, v) e M X (0, 1). Furthermore, v^(m, 1) = /" (pe^du =: H{u) and 
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Va(m,0) = 0. Application of the Laplace transform in u and since XH{X) = <p(/L — 1), the 
solution can be explicitly calculated as (cf. (6.3)-(6.5)) 

, sin(e(A)v) . , 1 sin(e(A)v) , 
sin(e(A)) A sin(e(A)) 

When taking derivatives of (6.29) in v, additional factors of eX are created; furthermore 
the multiplier has either the cosinus or sinus in the denominator. We have for v € (0, 1), 



sin(eAv) 



sin(eA) 



< 



cos(eAv) 



sin(eA) 



< 



C 

— ( 



,e|A|(v-l) 



< 



c 

en 



-e2 



ie|A|(v-l) 



(6.30) 



Indeed, for |eA| < 1, we have M = liSlijI < 



c 

\eX\ 



^ and lel^^K"-!)! > g-i. For 



(6.31) 



(6.32) 



\eX\ > 1, we have jU = ^ (and hence ^ = 1) |cos(eAv)| < 4el'^'^l^' and |sin(eA)| > ^e^'^^K 
This proves (6.30). Now, we multiply (6.29) by |Ap'^+';u'^, apply the norm on the 
line 9tA = 3£ — ^ and use (6.30). This yields the estimate and in particular uniqueness. 
This yields (6.26) for (note that 1 < |A| and 1 < /i for the considered values of A). This 
concludes the proof of (6.26). The estimate of (6.27) follows directly by multiphcation of 
(6.29) with |A| and evaluation at 9?(A + 1) (corresponding to one derivative in x). □ 

Lemma 6.6. Let ^ be the coordinate transform from Lemma 6.5. Then P : K'^ -^M.is a 
solution of (6.33) if and only ifp = Po^:K^R satisfies 

Aep = R{p,(p) inK, 
P = fx on diK, 

Py = on doK, 

where the operator R{p, (p) is given by (using the notation 7 := (1 + Yy)~^) 

R{P,(P) = -YxWxPy-YWxxPy-^YWxPxy + YYyVxPy 

+ fWxWxyPy + fWxPyy + ilfrVyPy + Hfif ' ^)Pyr 

Proof. We can write the inverse coordinate transform ^ ■.K'l'^Kas ^{x,y) = {x,y + 
Tj(jc,>')) for some 77 -.K^^K. In particular, y + 77 + + 77(1, y)) = j. By differ- 
entiating this equality ini andy, we get ^|fx + {\ + '^fy)v\x~^ and (1 + + 77-,-) = 1. This 
imphes 55 = 1> 3I = 0, 3^ = - JWx and j| = 7, in particular, P^ = 0-;^Py = Q, justifying 
the boundary condition in (6.33). Equation (6.32) follows from (6.33) together with 

Pxx = {dx-Wtdy){Px-WxPy) 

= Pxx- YxVxPy - WxxPy - 27 VxPxy + YYy VxPy + fVxVxyPy + fVxPyy, 
Pyy = Y{YPy)y = fPyy + YYyPy 

□ 

6.3. Estimates on the pressure. The main result of this section is: 
Proposition 6.7 (Shape dependence of p). Let k gN with k > I, £ G {0, 3(2fjr[)) '^nd 
(peXg /G-^I- Then there is a constant a < lsuchthatif\\^\\k-i/2 < (X, thenthere 
is a unique solution P -.K ^Rof 



AeP = 
P = fx 

Py = 



inKf, 
on diKf, 
on doKl>. 



(6.33) 
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Furthermore, p = P satisfies 

WpWyI < CWfWxl- (6-34) 

Let q be the solution of (6.1) with data f and with right hand side g = Q. Furthermore, let 
P, p and q be the corresponding solutions with boundary condition f instead of f and let 

w = p — q and w = p — q. Then 

lb-/5||j.| < C||(p-^||^,_,/2(ll/llx| + ll/llx|)+C||/-/llx|- (6-35) 

||w- wllj,, < C||9 - + +C||/-/llz|(ll9ll4-i/2 + ||9ll4-i/2). 

Before we address the proof of Proposition 6.7, we give an estimate of the nonlinear 
term R(p, (p) defined in (6.32): 

Proposition 6.8. LetkeN, k>l, e e (0, j^f^)- Let (p,^ eXg^^^^ with ||^||^t-i/2 < 

1 1 ^ 1 1 t- 1/2 < Q for some sufficiently small universal constant c^. Let p be the solution 

of (6.33) on K'^ with boundary data f, let p ^Yg be the corresponding solution on K'^ 
with boundary data f. Then 

||/?(P,(P)-/?(P,^)|U < C||/J-p||y.(||(p|| ,_V2 + ||^ll;,*-l/2) 

(6 36) 

+ C||<p-^ll;f|-l/2(||p||yt + ||p||y|), 

where R{p, (p) is defined in (6.32). 

Proof. The proof uses some ideas of the proof of the algebra property in Lemma 5.4. There 
are two differences: We also need to take a supremum in the angular variable v. The factors 
in R are controlled by different norms. Furthermore, the norm g controls a discrete set of 
homogeneous norms (cf. (2.38)). We will show that 

\\R{p,9)\\zl < C\\(p^-My, (6.37) 

the proof of (6.36) follows by an analogous argument, using the multilinear structure of R. 
We use the representation (6.32). Recall the definition (6.28) of the coordinate transform 
Xj/. We claim that (6.37) follows by iterative application of 

llp^llz* + ll(l-rkllz*< C||(p||^*-,/2||^||z*, (6.38) 

IMy^ + WC^PyMzl < C\\(p\\^,-,/2{\\w,\\^, + \\'-Wy\\^,), (6.39) 

where p is either one of the terms e i/'x or i/a^, and where we recall that 7 = ( 1 + i/Ay ) ^ ' . 

Assuming that (6.38)-(6.39) hold, the proof of (6.37) is easy: Recall that with the defi- 
nitions of the norms (2.34), (2.38) and by (6.34) we have 

b^llz| + ll^/'yllz| + ll/'«|lz| + ll^/'xy|lz| + ll(i)'Pyjllz| < C\\p\\yk. (6.40) 
We show how the estimate of (6.37) proceeds for the term YxVxPy (the first term on the 
right hand side of (6.32)): In view of Yx = —Y^Wxy, we have 

, , (6.38) (6.39),(6,40) 

WYxVxPyWzl = \\r{eVx)Vxy{lPy)\\zl < <^ll V^;^(b)') llz| < <^ll<Pll4-i/2||/'lk^ 

Indeed, it can be easily checked that the estimate of the other terms in (6.32) proceeds 
analogously. In order to conclude the proof of the Proposition, it hence remains to give the 
argument for (6.38)-(6.39). 
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Proof of (6.38); In view of the Taylor expansion y— 1 = y/y — 1//^ + . . ., the estimate for 
the term (1 — 'Y)g in (6.38) hence follows by the corresponding estimate for the term \\lyg 
together with the estimate (6.26), ||<j!)|| *-i/2 < and for q sufficiently small. In order to 
see (6.38), it hence remains to show 



(6.41) 



Let Vxff and Vg be the Taylor polynomials at {x,y) = (0, 0) of and g (of order 3A: — 1 and 
Tik — 4, respectively). We decompose i/A = i/Ai + y/o and ^ = + where i/^i := V\j,tl, and 
gi :='Pg^ and where the radial cut-off function = e C"([0,°o)) satisfies C = 1 

in [0, g] and ^ = in [| , «>] . The terms po and pi are defined correspondingly. Clearly, it is 
enough to show the corresponding estimate to (6.41) for the products pigi, Pigo, Pogi and 
pogo- In the following, we will give the estimate for pogo', the estimate of the other terms 
(related to finite dimensional Taylor expansion) proceeds analogously as in the proof of 
the algebra property of X| in Lemma 5.4. Furthermore, for simplicity of notation, we give 
the proof with A replaced by |A| in (2.37), i.e. when only radial derivatives appear . The 
argument in the case of angular derivatives ^dy proceeds by distributing the derivatives on 
the two factors using Leibniz' rule. We will use the notation ^g'^' =e2^l'*'l(^~^), in particular 

\X\ lyL— nl It?I 

by the triangle inequality we have Kir ' < K'^ "TT^". 

Analogously to (5.6), we decompose the functions ij/o and go into their low and high 
frequency part, i.e. Yo = V+ + V- (with p± defined correspondingly) and = ^+ +^-- 
In particular, as in (5.7), we have 



|Anf;i + (i)*lir| < Cm* 



|A|*|p;i + (i)*|pr| < Cm^IPoI. (6.42) 



We need to estimate the terms g+p+, g+p-, g-P+ and g-p-. We show the estimate for 
the high frequency /high frequency product g-p-. Indeed, the estimate for the other two 
terms proceeds similarly (see also the algebra proof in [26]). 
In view of (2.37) and (6.42), we need to show for all £ e [hk], 



\P-g- 



C\\ sup I / K, 
v€(0,l) •^9?i7=r 

lA\g\\zp 



< c\W\\^,-, 



p-{^-n)g-{n)(^{n)\\y[m=u-i) 

(6.43) 



for any 7 e M of our choice such the above integral is defined. Let JC G [0, 1) be the 
smallest number such that A^:=3^ — 2+K"eNo. We show the corresponding stronger 
estimate to (6.43) where the A^*^^^ is replaced by . The estimate is stronger since |A | > 5 
in the line of integration in (6.43). The advantage is that the binomial formula = 
1^=0 CVXA-T/)-'?]'^--' withAr = 3^-2eNo canbeappUed. We have 



L 



sup 

9?A=3^-5v6(0,I) 



-L 



9?r]=7 



KpX 



P-{X-T])g-{T])(B{T]) 



cB{X) 



sup 

9?A=3^-5vG(0,I) 



i=0 ■'^ 



U-2+K-i 



P-{X-ri)ri'g-{ri)cB{ri) 



cB{X). 



Note that the above inner integrand is analytic in 77 and hence does not depend on the value 
of 7 as long as the integral is well-defined; hence we may choose 7 freely as a function of 
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i. In turn, since furthermore k'^^^ < K^g ^^k}^ \ it is enough to estimate terms of the form 



/ sup / /^i^-''l|A-T,r-2-,|^(A-T,)|4''l|T/|i|r(r?)I^S(Tj) 



2 

d3(A) 
(6.44) 



for £ e [5, A:], all integers i € [0,3£-2 + Jc] and with our choice of ye M. We next apply 
the following variant of (5.10) which says that for all 5 > j, we have 

II sup |F*G|||i2(5R;^=^) < Qll sup il + \Mh\F\\y(^^^^pJ sup |G||L2(9,;t=&) 

V6(0,l) ^ V6(0,l) ^ ' V6(0,l) ^ 

(6.45) 

if J3i + J32 = /3 and as long as all integrals are well-defined. We introduce the short notation 

ll^lkA=j3 = II sup^g(o,i)^1^'|^IIIl2(9{A=;8)- In view of (6.44) and (6.45), for the proof of 
(6.43) it suffices to show for all ^ e [ j , A;] and for all integers i e [0, 3£ - 2 + fc], 

ll(l + |Ap'+'^-'-'+^)|p=|||m=ftll|Ar^||m=ft < ll<Pll4-i/2|kllz|, (6.46) 

where we can arbitrarily choose 5 > j and /3i , ^2 with J3i + /Ji = 3^ — 5 . With the notation 
5 = 5 + (jc — 1) < 5, it is equivalentiy enough to show . 

||(l + |A|3^-l-+^)|pr|||95;t=ft|||Ar-fr||cR,=ft < ||<p||^.-V2|k||z|, (6.47) 

where we can arbitrarily choose 5 > 5 and J3i,j32 with j3i + j32 = 3£ — j. Both j3i, P2 as 
well as 5 are allowed to depend on i and i. In fact, in the sequel we will always choose 

Recall that either po = exj/ox or po = Voy and hence either |po(A)| = \X\i/o{X + I)| or 
|p()(A)| = \dv\j/o{X + I)|. By (6.26) and by the same argument as in the proof of Lemma 
5.2, we hence have 

||(I + |AP^+i/+')IPo|||<r;l=3, 1 < C||(p|| V£G[0,A-i]. (6.48) 

By (6.48) and (6.42), we have 

||(l + (ir'|Ap^+^)|^|||^,^3,_, < C||(p||^,-v2 V£e [0,A-i]. (6.49) 
Furthermore, in view of (2.37) and by (6.42), we also have 

\\CsY\M''-'\g^\hx=M-l <c\\g\y^, y£G[l,k]. (6.50) 

We prove (6.47) for the three extreme cases, i.e. the corners of the triangle in (ij) where 
i and / e [0, 3i? — 2]}, the estimate of the other terms follows by 'interpolation' of 

these estimates: 

(al) £=j,i = (and hence j3 = 3^ - | = -|): With the choice j3i = 3 • i - i = 0, 
j32 = 3- | — 5 = — 5 and 5 = j, we need to estimate 

||(I + |A|i)|pr|b,=oll(^)^nikx._3 < C||9||^|-i/2||g||4. 

(bl) £^k,i = 0: With the choice J3i = 3(A;- i) - ^ = 3A;- 2, J32 = 3 • ? - ^ = -| and 



5 = 5 , we need to estimate 



2 J 2 ~ -"^ ^>P2 — -'3 2 ~ 2 

1 

||(I + |Ap'=-2(i)*-t)|pr|||35,=3,_2 ||(i)t|^|||3,,_3 < C||<p|| ._V2|k|U 



1 , r ,;,_2 , ,„ ,,,1,2, 
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(cl) i = k,i = 3k-2: With the choice j3i = 3 • ^ - ^ = 0, jSi = 3;t- ^ and 5 = i, we need 
to estimate 

||(l + |A|i)|pr|||5R;t=oll|Ap*-2(i)*inilm=3;t-J < CM^,^-Mzi- 

Now, using (6.49) and (6.50), it can be easily checked that the above estimates (al), (bl) 
and (cl) hold true for k> I. This concludes the proof of (6.41) and hence of (6.38). 

Proof of (6.39).- The proof of (6.39) proceeds similarly to the proof of (6.38). As before, 
we show the estimate for the crucial high frequency/high frequency case p-^w-. With the 
same arguments as before, we need to show (correspondingly to (6.47)): 

(|)*ll|Ap'-'-'+'^|p^||kA=ftll(l + |Ar'+^)|w^|||m=fe 

< C\\(p\\^,.y2{\\w,\\^i + \\lwy\\^i) =: n, (6.51) 

(we have put the 5 on the second factor). In the above inequality, we can arbitrarily choose 
j3i , j32, 5 as long as j3i + j32 = 3^ — j and 5>\. Both /3i, ^ as well as 5 may depend on 
I, i. We note that by (2.37), we have 

\\{\Y\M''-'\^-\\Wx=u-l <C{\\w,\\zi + \\l^y\\zi)^ V£e[i,fc]. (6.52) 

We prove (6.5 1) in the case when the maximum number of derivatives fall onto p, i.e. i = k 
and i = (in particular, K = where we recall that k is defined as the smallest nonnegative 
integer such that 3£- 2 e N). With the choice j3i = 3{k- j) - 5 = 3^:-2, /Jj = 3 • 1 - | = ^ 
and 5 = 1, we hence need to estimate 

ll(i)*|Ap*-i|pr|||c«,=3*_2ll(i + |A|)(i)i||r||i3,;^^. < cn. 

This estimate holds true as can easily be checked using (6.49) and (6.52). The estimate of 
the terms P-,xW-, P-^W- and p-,xW- proceeds similarly. This concludes the estimate of 
(6.39) and hence of the proposition. □ 

Proof of Proposition 6.7. We first show the existence of a solution P of (6.33) on K^. By 
Lemma 6.6, for this we need to find a solution p of (6.31). We will solve (6.31) using an 
iterative argument. We set po := and iteratively define to be the solution of (6.1) 
with right hand side g = R{pi-i , (p) and boundary data fx. By (6.2) and (6.36), we get 

(6.2) (6.36) 

WPi+i-PiWri ^ C\\RiPi,(p)-R{Pi-i,(p)\\zl < C\\(p\\^,-i/2\\pi+i-pi\\Yk 

< Ca\\pi+l - piWyk < \\\Pi+l - PiWyk 

for a sufficiently small, {pij/gN is a Cauchy sequence and converges to a solution p of 
(6.3 1). By (6.2), p satisfies (6.34). The estimates (6.35) now follow from the representation 
(6.31) of the solution together with the estimates (6.36) and (6.37). □ 

We also have the nonhnear version of the trace estimate in Lemma 6.4: 

Lemma 6.9. Suppose that the assumptions of Proposition 6.7 are satisfied (in particular, 
k>\). Then with the notation of Proposition 6. 7, we have 

Proof. Indeed, (6.53) follows by appUcation of (6.17) and since p satisfies (6.31). □ 
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Note that the corresponding 'bilinear' estimates for w, p — p and w — w also hold (cor- 
respondingly as in Proposition 6.7): The estimates for the case of x derivative are 

l|5.K)|r|l4-i < C||(p||^,_v2||/||4, 

I|5.(p.)|r-5.fe)|rll4-i < C||(p-^||^|-V2(ll/llz| + ll/llx|)+C||/-/l4, (6.54) 

||5^K)|r-5^(w;c)|rll4-i < C\\(p-(p\\^,.,,2{\\f\\^k + \\f\\xi) 

+ C|l/-/ll4(ll'Pll;,^*-V2 + ||9||^*-l/2). 

These estimates follow analogously lo (6.53) but by using the 'bilinear' estimates (6.36) 
for R. The corresponding estimates to (6.54) with dx{wx)\Y and dx{px)\T replaced by 
{\y'dx{wy)\Y and (^)^^jc(wy)|r. respectively, also hold. 

6.4. Estimates for the profile. In this section, we prove Proposition 4.2. We show that 

\mf)\\Txi+\\m)\yix2) < cmi^,^,. (6.55) 

The proof of (4.3) then follows by a straightforward extension of this estimate (indeed the 
estimates (6.36) show that the main part dxBgdx of the operator Af^ can be estimated like a 
biUnear operator). We recall, 

Neif) = ^Vrf^^^n^) +^Xf -f) +Aef, 

(l + e2)iV l+e^ Jix=o V(l+e2(i+/)2)i; 
cf. (2.18). The two main estimates in favor of (6.55) are: 

\\fxx\x=ofx\\TXl + \\fxx\x=ofx\\L2^Xl) ^ ^11/11^4+1' (6.56) 
\\dxB{dxf+Aef\\j.x, + \\dxB{dxf+Aef\\^^, < C||/||2 (6.57) 

Two "nonlinear" corrections have been neglected in the estimates (6.56) and (6.57): One 
correction is related to the —3e^ f^\^^^Q term, the other correction is related to the term 
fx{{l + s^{l + f)^) — I). Indeed, these corrections are easily controlled as lower order 
terms for sufficiently small /; the estimate of these terms is left to the reader. In the 
following, we will give the proof of (6.56) and (6.57). 
Proof of (6.56): In view of the definition (2.30), (6.56) follows from 

[^/(/«|.=o/x)]i2(;,.) < Cll/ll^^,^, , forO<i + j< k. (6.58) 

In order to see (6.58), we note that by (5.2), we have 

\\fxx\x=0fx\\x0 = \\fxx\x=0fx\y < \\fx\\L2\\fxx\\L- ^ C||/||^V2 , (6.59) 

IIA.|.=o/.||;,; < IIMI;^; Wfxxh- ? 11/11^,+, for l<j<k. (6.60) 

We take the square of (6.59) and integrate the equation in time. We then apply Holders 
inequality with on the || • term and LT on the \\ ■ \\ 1/2 term. In view of (5.12), this 

yields (6.58) for i — j — 0. The estimate for / > or /' > follows from (6.60): We take the 
square of (6.60) and integrate in time. We then apply Holder's inequality with on the 
II • II j+i term and L°° on the || • || j term. In view of (5.12), this yields (6.58) for i = and 

Xe Xg 

1 < j < k. It remains to prove (6.58) for / > 0: In this case we note that the time derivatives 
on the left hand side of (6.58) are distributed according to Leibniz' rule on the two factors. 
The loss of regularity due to the time derivation is compensated by the fact that we only 
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need to estimate the smaller norm j = £ — i. In view of the definition of || • this 
yields (6.58) for all < ? < ^ and 7 > 1 . This concludes the proof of (6.56). 
Proof of (6.57); With the same argument as before, the proof of (6.57) can be reduced to 
the following estimate which does not involve time: 



[d,BidJ+Aef]^j < C\\f\\ 



where we may assume that ||/|| k+1/2 < 0! < 1 (this corresponds to the estimate ||/||^, k-1/2. 

a < 1 for the corresponding time-space estimate). Let 7 be defined in Lemma 6.5 and 
let p, q be defined as in Proposition 6.7. By (2.17), (2.15) and (2.8), we have 

(l + e2)-ia,[^,-(i)2^3,],r, 

+/) {p, - YY^py) - {IfrPy] |r- 



l/ll 



-7+1 



for 0<j<k, 



(6.61) 



< 



Aef, 

-dMfx = (1 

With the notation w = p — q and Cg 

Aef + dXfx = Ced^ 
We first note that for aU 7 > 0, we have 

IIK|r).||;,7 + ll(i)'(w: 



(1 



< 1 , we hence get 



Wx-fPx + ii+f)7VxPy + {^)H7-i)P: 



(6.54) 

< 



^7+1/2 1 



1,7+1 . 



(6.62) 



(6.63) 



For the other terms on the right hand side of (6.62) and for y > 1 we additionally use the 
algebra property (5.8). The estimate of the term f{Px\r)x proceeds as follows 



(5.8) 



I|(/P.|r)x||;,7 < ||/.||;,7||p.|rll 



(5.3)^(6.53) 

||/||;,7||(;..|r).|l;,7 < C|l/|l 



,7+1/2 



l/ll 



-7+1 7 



where we also used ||f;c|r||L2 < C||p||y^i < C||/||^^i which follows from (5.11) and (6.34). 
We also note that in view of (6.27) we have ||ev^dr|L«:+i/2 + ||VA)|r|L*+i/2 < Ca < 1 if 

'^e X^ 

the constant a in the assumption of Theorem 3.1 is chosen sufficiently small. In view of 
the Taylor expansion 1 — 7 = Xj/y^ Vy + ■■ ■ ^nd by (6.27), this implies also |1 'J^r ~ 1 11^7 ^ 
II Vyirllx^ — ^- P^'^ 7^1' the estimate of the remaining terms in (6.62) then follows using 
these estimates together with (5.8), (6.27) and (6.53). 

It remains to give the estimate for j = 0, where the algebra property (5.8) does not 
apply: We show the estimate for the nonlinear term {fPx\r)x, the estimate for the other 
terms proceeds analogously. Indeed, we have 



ll(M|r).|lL2 < ll/lkHI(/'.|r).|lL2 + ll/.|lL2b.|r||L- 



(6.53) 
< 



C||/||^l/2| 



\\X'' 



where we also used ||/7f|i ||^2 < C||/5l|j,i < C]|/||-yi which follows from (5.11) and (6.34). 
The estimate of the other terms proceeds analogously. This concludes the proof of (6.61) 
and hence of the lemma. 



7. Localization argument 

We prove Propositions 4.4-4.5 thus concluding the proof of Theorem 3.3 in Section 4. 
We will use the notation used in the proof of this theorem. We first note that the derivative 
5£e of £e at w is given by 

= L^f + K^f, (7.1) 
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cf. (2.9), where the operators Lg (leading order) and Kg (remainder) are given by 
and 

K^f = (Z)"' - I) ft +lrf+d,{V"fw) (73) 
- [{{l-x)s-+xs+)f]^- [{{l-x)a^+xa^)w]^ (7.4) 
_ -^H,. 3w,{2e^w,f) 
- 2(l + e2(w + /*)2)i 

Here, we have introduced the following notation: defined as in (2.10) with / 

replaced by w; is defined as in (2.4) with s±{t) replaced by s^{t). Furthermore V^f = 
5D^{f) and (y±^wf = 5s±^w{f) are the linearizations of D^, a±^w Furthermore, 
denotes the shape derivative of (derivative in w direction). A dot on the top of a symbol 
denotes the time derivative. 

The proof is based on two small parameters 5 , T > 0. The parameter 5 is used to locaUze 
the estimates near the boundary. Note that in the interior where C > h > c > 0, the operator 
is uniformly parabolic. The parameter T signifies the time integral where the solution is 
defined. In the course of the proof, we will choose 5 and T = t(5), in this order, to be 
sufficiently small. In this section, we write c, C for all constants which depend only on 
/in, A:, but neither depend on 3 nor T. 

For 5 > 0, we define a covering of E = (0, 1) by setting E^g := (0,25), £25 '■= (5, 1 — 
5), £35 : = (1 — 25. 1). Correspondingly, let 2,5^ := (0, t) x Ejg. We choose a smooth parti- 
tion of unity { V!5}i=i,2,3, subordinate to this covering with xj/jg e C~(£,-5, [0, 1]), J^i ViS = ^ 
on E and also ||5-'Vr5 ^ Cd~j for all j > 0. We also define a second set of cut-off func- 
tions {V>i5},=i,2,3: Let £15 := (0,45), £35 := (35, 1 - 35). £35 := (1 -45,1). We choose 
V15 € C~(£',-5, [0, 1]) with xjfjs = 1 in Ejg and in particular, ^figij/ig = y/jg. Finally, we also 
consider the cut-off function T^ig, supported on an even larger set such that ^isWiS = Wi5- 
The support of y/ig = 1 is e.g. included in (0,55), the fimctions V'25)V'35 are defined 
analogously. 

Since ll/inll-^s/z^^^ <Cand j-^i^g^j < C, it follows that /m and w are Holder continu- 
ous in time and space. Therefore, by the boundary condition |/zin,x(0)| = |/zm,x(l)| = 1 and 
recalling that hin = x(l — x) + f\a{x')dx' , there is 5o > such that for all 5 < 5o we have 
\dxhQ\ G (2,2) in E^gdE-^g and ho S (c,C) in £25- Similarly, there is 5o and To such that 
the corresponding estimates hold for h^{T) = x{l —x) + jQw{T,x')dx' for all 5 < 5o and 
every fixed time T < Tq. In the sequel, we will always assimie < 5,T < min{5o,To,0.1}. 

7.1. Proof of Proposition 4.4. 

Proof of Proposition 4.4. In view of the extension Lenrnia 4.3, we only need to consider 
the case of zero initial data so that in the following we may assume fi„ = 0. We begin with 
the proof of maximal regularity estimate (4.13). That is we assume that / satisfies 

5£e/ = gin<2T and f = on(0,T)x^£ (7.5) 

and with initial data / = and we will show 



+ 



5B-(/, 



(7.6) 
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With the partition of unity xj/jg, i ~ 1,2,3 and by the triangle inequality, we have 

ll/llrx|+'(eT) - \\^lsf\\TX^+i^Q^) + \\W2Sf\\TX^+'{Q,) + \\^iSf\\TX^+i(Q^y (7.7) 

We begin with the estimate for y/isf (related to the left boundary of the domain). The idea 
is to use that on Qist (that is near the left boundary), dC (and also Lg) are approximated 
by Lg, where L^f := ft + Ag/ and where Ag is defined in (2. 17). We first claim that 

WWisfWrxt^Q.) ^ CsWMVfMTxIiQ,) +^5 11/11 rx|(G.) + ^WfWTxt'iQ.y ^^-^^ 
where we recall V''i5V'^i5 = ¥i5 ^^d hence (1 — V'^i5)v^i5 = 0. In order to see (7.8), we first 
note that by (4.2), we get \\Wi5f\\TxtU.Qz) - *^II^'^V^i5-^II7'x|((o.t)xa:) (*^ ^-^-term is esti- 
mated by the other term since the support of the right hand side is bounded). Furthermore, 

LeiWlsf) = ¥lsLei¥lsf) + {^-¥ls)LEiVfisf) 

= WisMWisf) + Visft + (1 - Vi^i5)Ag(v^/). 
Furthermore, by (2.17) and since xj/ig depends only on x, by a short calculation we obtain 

M¥l5f)=dx{-p'x+^Py), 

(1 - VAi5)Ag(VAi5/) = dx{-p'^ + ^p'y) - d^fisp' - 2dxfi5p'^ + ^dxWl5Py, 
where p' -.K ^M. and p' := (1 — ^H5)p' the solutions of 

A,p' = inK, ( Aep' = -Id^Yisp'x-d^VisP' ^K, 

p'^iWisDx onr = diK, I p' = ondiK, (7.9) 

Py = on ^0^- [ = on doK. 

In the following, we use the following ReUich-type estimate for lower order terms which 
can be obtained by standard interpolation: For any given 5 > 0, we have 

II/X.II4-1 < 5 ll/ll^i + Q 11/11^,-1. (7.10) 

Using this inequality, estimate (7.8) follows by application of Proposition 6.1. 
Let [Lg , v^is] := ^WiS ^ WisL" denote the commutator of Lg and i/Zig, then 

\\¥ldLeiWlsf)\\TX^^(Q,) < Wisii^e - SC^^){fWls)\\TX^(Q,) + \\W\5S^{WlS.f)\\TX^(Q,) 

< WisiLe - 5/:e)(M5)llrx|(eo + WidiSC'^, WMtx^q,) + WWissWrxliQ,)' 

(7.11) 

where we have used WidVi8 — V\5- Th^ estimate of the first and second term on the right- 
hand side of the above estimate is given in Lemmas 7.1 and 7.2. Choosing 5 sufficiently 
small, we hence obtain 

WWldfWrxt'iQr) - ^s\\s\\TXl{Q,)+Cs\\f\\TXl(Q,) + l MtX^+^Q,)' 

The third term on the right hand side of (7.7) can be estimated analogously. For the middle 
term (corresponding to the interior of the domain), we note that in the interior, our weighted 
Sobolev norms are equivalent to standard Sobolev norms (with equivalence depending on 
5) and an analogous estimate to the one above can be achieved for the middle term using 
standard parabolic estimates, see also [16, 29]. Altogether, these estimates yield 

(7.7) 1 

11/11 rxi+Hex) - '^slkllrzKe^) +C5ll/llra|(e^) + 2"-^"^^e+ne0" 
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Then using that /j;=o = 0' we deduce that Il/ll7'x*(e^) — ^'^^''^ll/llrA'*+'(2x) ^® choose 
T such that QCt^/^ < g. Hence, 

which yields (7.6) by absorbing on the left hand side. 

The existence part is similar to the proof of Lemma 3.4 of [20]. Indeed, the argument 
used there can be generalized since only very little of the particular structure is used: The 
main ingredient in this argument is existence and maximal regularity for the linearized 
operator at the boundary. The second ingredient is the fact that existence of a solution to- 
gether with estimates in the interior follows by standard parabolic theory. We have already 
proved these properties for our operator. Finally, the argument also requires that the oper- 
ator can be locaUzed in the sense that the long-range interaction of the solution operator 
only yields a lower order contribution. Indeed, we have used this idea already in the proof 
of (7.8). □ 

In the following we give the estimate for the right-hand side of (7.11). We use the 
notations and assumptions of the proof of Proposition 4.4. 

Lemma 7.1 (Estimate of difference). For < T < 1, we have 

Wf.siLe - W)llr4(e.) < Qll/llrxKe.) +C'5||/|lrx*+i(2.)- ^'^'^^^ 

Proof. Let p' and P be the solutions of 

Aep'=0 inK, ( AgP = in £2"', 

p'^{Visf)x ondiK, I P={Visf)x on^iH"', (7.13) 

p'y = ondoK, [ Py = ondoSr^, 

where n*" = {{x,y) |0<x<l,andO<y< and diST = graph/i"" and /z"' = x(l - 

x) + fQw{x')dx'. The existence of a solution p' follows from Proposition 6.1, the existence 
of a solution P can be shown with similar arguments as in the proof of Proposition 6.7; 
they will not be detailed here. The reason to introduce these two functions is that we 
have Le-L^=Ae- A^. Furthermore, Ae{\l/uf) = d,{-p', + ^p'y) and A^(v^i5/) = 
dx{—h^Px + jiPy)- Note that A^ and A^' are not defined on the same interval in x. To 
compare them, we therefore use the cut-off function y/jg. Indeed, "fi^P can be seen as a 
function in the domain since the domains and £1^ coincide on the 

support of ^isP (recall that "^\s^\s = Vfis)- More precisely, ^figP solves 

AeiflsP) ^IdxYisdxP + d^flsP in/:V>'is("'-f), 
iWisP) = {Wi5f)x ondiKVisi^-i), (7.14) 

{yf^gP)y = on doK'^isi^-i) . 

We use Lemma 6.5 to construct a coordinate transform W = id) from K to K^isi^-i)_ 
Hence, arguing as in the proof of Proposition 6.7, we deduce that p = Po^ solves 

AeiVisp) = R{p,V)+}<^i{P,V) in^, 
TT = ondiK, (7.15) 

Tly =0 on doK. 

where R{p, is given in (6.32) and where /Ci {p, yf) is a lower order term involving at most 
one derivative of P and such that supp/Ci C supp(9;c\/Ai5. Arguing as in the proof of Propo- 
sition 6.7, we deduce that ^fisp S and hence ^fisP S Fg (O"'). Similarly, it can also be 
shown that {yf^s + W35)P (^"'). One can then start a bootstrap argument and prove that 
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P e y^Hii'"), WPWy^ia^) < C(w)||v/i5/|lx| and \\P\\yM^n„^ < C(w) ||v/i5/ilx|+i- Note 
that yffisP' satisfies the same system as (7.14) when replacing P by p' and ^"?^i«("'~2) by K. 
Taking the difference, n = p — p^we hence obtain 

Aeifis^t) = R{p,V)+^2{p,p',V) in^, 
TT = on diK, 

TCy = on ^0^- 

where K2{p,p' , ^r) is a lower order term with at most one derivative of por p' . We have 

V^15^nri5/) = ( - K{Wl5P)x + ^{WlsP)y + h:P{Wls)x), - dM-h^P, + ^Py), 

WisMWisf) = ( - i¥isP')x + ^ ("KisPOy +p'{Vi5)x)^ - dxWisi-p'x + ^P'y)- 
Hence and since — 1 1 < C5 on supp i/z^ig, we obtain 

= {-{¥l5^)x+^{Vl5^)y + ^{¥l5)x)^-dxWl5{-^x+^7Cy) + K:3, (7.16) 

where the remainder term /C3 satisfies ll^3|lrx|(eT) — ^^ll/llrz*+'(eT)- Notice that the 
second term on the right-hand side on (7.16) consists of lower order term. Hence, we only 
need to estimate the first term. We have 

||/C2||2*+i <C5||/||;,| + 5||/||^.+, and \\R{p,w)\\^,^^ <C5\\f\\^,^i. 

Hence, ||7r||y<:+i < CgH/Hj^* +C5||/||^t+i and (7.12) follows easily. In a sense we proved 
that the difference between A* and comes from either terms which are small or terms 
which are more regular. This shows the estimate of — L", the estimate of follows 
similarly also using (7.10). This concludes the proof of the Lemma. □ 

Lemma 7.2 (Estimate of commutator). For < T < 1, we have 

||v/i5[5££,v/i6]/ll7'x|(eo ^ '^5ll/llrx|(ex)+'^^ll/llrx|+i(er)- 
Proof. Let P and Q be the solutions of 

A>;P = inil"', ( AeQ = inil"', 

p = {vi5f)x on^l^^ \ Q^fx onal^2^ (7.17) 

Py = ondiQ'^, [ Qy^O on^o^"'- 

Arguing as above, we can prove that Q,P € Y^+^{Q.'*') and that ||2||y*+i(Qw)> ||^||yt+i(£2>v) < 
C(w)||/|Lt+,- for i = 0, 1. Moreover, 

¥l5[Le,Vfl5]f = Vfl5[{-WxPx+^Py)x-Vfl5{-WxQx+^Qy)x] 

= Wl5[{-MP-¥l5Q)x+^{P-WlsQ)y)x]+}C, (7.18) 
where the terms of lower order are collected in /C. Note that P—^isQ solves 

Ae{P - W15Q) = -2dx\l/i5dxQ - d^Xj/is Q in Q.^ , 

P-Yi5Q = ¥i5.xf ondiQ", (7.19) 

{P-¥i5Q)y = ondoQ.^. 

The right-hand side of (7.19) is more regular and allow us to estimate P— ij/idQ by 
11^ ~ VisGllyife ^ ^ll/llx* + To • Moreover, the operator K on the right hand side 
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of (7.18) only involve one derivative of Q and hence can be estimated similarly. This 
concludes the proof of the Lemma. □ 

7.2. Proof of Proposition 4.5. A localization of the estimate in Lenmia 5.6 yields the 
following estimate: For all A;, y e No, T > and /o e (Qt) we have 

l|5//0|lc0(4+l/2(£„ < C(||^/+V0|li2(;,|(£)) + ||a//0||^,(^.+.(^))). (7.20) 

In view of Lemma 5.6, it remains to show boundedness and continuous differentiability of 
Cg, defined in (2.9). The highest order term of £e is given by the operator dxB{, see 2.8. 
Boundedness of this operator follows from a localization of the estimates in Proposition 
6.7 and Lemma 6.4: For p be defined by (2.8), we have 

WdxBiMxIiE) < C(||a.(;>.)|r||x|(£) + (^)'l|5.(P.)|r|lx|(£)) 

where § is a lower order term. Indeed, near the left boundary of E, g is the same term as 
in the right hand side of (7. 15). In the center and near the right boundary of E, g is defined 
analogously. By a localization of (6.37) and since — hence obtain 

Furthermore, also using that — ^snq 

¥±{f+n\\xi < ciMiik- (11/114 + iirii^) < cii/ii^: + iirii^.) < cii/ii^.+. . 

It remains to show boundedness and continuity of the first derivative 5£g : In view of 
(7.1), the estimate of 5£g leads to analogous terms as the one for the estimate of £e- The 
estimate hence follows analogously and is not more difficult than the estimate of Ce itself. 
Similarly, continuity of the first derivative follows by a bound on the second derivative. As 
before, in view of the structure of the operator it is clear that these terms do not impose 
any other difficulties than the one's already when estimating the operator Cg itself. Hence, 
a bound on the second derivative can be obtained similarly to the calculations before. 
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